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Abstract 

We report on the results of a spectroscopic survey of 130 comets that was conducted at McDonald 
observatory from 1980 through 2008. Some of the comets were observed on only one night, while 
others were observed repeatedly. For 20 of these comets, no molecules were detected. For the 
remaining 110 comets, some emission from CN, OH, NH, C3, C2, CH, and NH2 molecules were 
observed on at least one occasion. We converted the observed molecular column densities to 
production rates using a Haser (1957) model. We defined a restricted data set of comets that had 
at least 3 nights of observations. The restricted data set consists of 59 comets. We used ratios 
of production rates to study the trends in the data. We find two classes of comets: typical and 
carbon-chain depleted comets. Using a very strict definition of depleted comets, requiring C2 and 
C3 to both be depleted, we find 9% of our restricted data set comets to be depleted. Using a more 
relaxed definition that requires only C2 to be below a threshold (similar to other researchers), we 
find 25% of the comets are depleted. Two-thirds of the depleted comets are Jupiter Family comets, 
while one-third are Long Period comets. 37% of the Jupiter Family comets are depleted, while 
18.5% of the Long Period comets are depleted. We compare our results with other studies and find 
good agreement. 

Keywords: Comets, composition. Spectroscopy 

1. Introduction marked difference being the amount of contin- 

uum present. This simple observation leads to 
Since the time that the first spectra of comets the question of whether all comets share the 
were obtained, it has been remarked that all same composition or if there are comets with 
cometary spectra are very similar, with the most fundamentally different compositions. If differ- 
ences are seen, then an additional question of 

the origin of the differences would be raised: are 
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questions are not just curiosities. We believe 
that comets are leftovers from the formation of 
the Solar System, that they have undergone only 
small changes since their time of formation, and 
that they represent the building blocks of the 
jovian planets. Thus, understanding their chem- 
istry is important for our understanding of the 
conditions at the start of the formation of the 
Solar System. 

With these questions in mind, we began a pro- 
gram in 1980 to obtain spectra of comets us- 
ing low resolution spectrographs at McDonald 
Observatory. This program continued through 
2008. We obtained observations of 132 differ- 
ent comets, some at more than one apparition. 
This paper reports on the observations of 130 of 
those comets. 29P/Schwassmann-Wachmann 1 
is left out because its almost-circular, 5 AU orbit 
and random outbursts make it an unusual comet 
that has been described in several previous pa- 
pers (Cochran et al. 1980, 1982, 1991a, 1991b; 
Cochran and White 1993). D/Shoemaker-Levy 
9 broke into many pieces prior to impacting 
Jupiter in 1994. Observations of the pieces had 
to be obtained in a different manner from the 
other comets and these are detailed in Cochran 
et al. (1994). 

Some of the comets described in this present 
paper were already discussed in prior papers (e.g. 
Cochran et al. 1987, 1989, 1992, 2009; Cochran 
and Barker 1987). They arc included again for 
completeness so that all the observations can be 
found in the same place. In addition, all obser- 
vations were analyzed in a consistent manner in 
this paper. Thus, this paper supersedes its pre- 
decessors. 

This is not the only survey of its kind to 
be carried out. A'Hearn et al. (1995, here- 
after AH95) reported on observations of pho- 
tometry of 85 comets. Fink (2009, hereafter 



F09) presented spectroscopic observations of 92 
comets. Langlund-Shula and Smith (2011, here- 
after LSll) observed 26 comets. We will com- 
pare our results with these studies near the end 
of this paper. 

2. Observations 

The observations described in this paper were 
obtained with three different instruments on two 
telescopes at McDonald Observatory. The vast 
majority of the observations were obtained with 
two spectrographs at the Cassegrain focus of the 
2.7m Harlan J. Smith telescope. A few observa- 
tions were obtained at the Cassegrain focus of 
the 2.1m Otto Struve telescope. 

The earliest observations were obtained with 
the Intensified Dissector Scanner (IDS) spectro- 
graph on the 2.7m telescope. This instrument 
had two 4x4 arcsec apertures, separated by 52 
arcsec along an east/west line, that were imaged 
on a photocathode. Spectra were obtained by 
rapidly scanning the complete wavelength range, 
which was typically 3200 - 6OOOA, at llA reso- 
lution. A complete spectral scan took 50 sec 
so all wavelengths were essentially observed si- 
mT.iltancoT.isly. The comet optocenter was alter- 
nately imaged through the "A" slit with the "B" 
slit imaging 52 arcsec away in the coma. Then 
the telescope was nodded to image the comet in 
B and coma in A. The seqTTcnce generally went 
ABBA. Sometimes, the telescope was set so that 
the comet optocenter was between the slits or 
at some other known position and spectra were 
obtained without nodding the telescope. Scans 
were repeated until sufficient signal/noise was 
obtained. 

The photocathode made the instrument ex- 
tremely sensitive. Thus, we were able to obtain 
spectra of comets as faint as V=19.5. However, 
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it also made it impossible to observe very bright 
objects. Thus, our standard stars were limited 
to those that were fainter than V=10 and were 
taken from the list of Stone (1977). In addition, 
we were not able to obtain observations of solar 
analog stars with this instrument (none are faint 
enough) and had to use fluxes from the Arvesen 
et al. (1969) atlas convolved with the slit func- 
tion of the instrument. The Arvesen et al. atlas 
is a whole-disk integrated high-resolution atlas 
suitable for our purposes. It is not the only avail- 
able atlas and there are known differences in the 
flux calibration between the Arvesen et al. atlas 
and other atlases, especially below 4000 A. Thus, 
the Arvesen et al. atlas might not be a perfect 
match for the Sun. However, because we color 
weight the solar spectrum (see discussion below) 
before we use it, the most critical feature for our 
purposes is that the lines are correct in wave- 
length and strength. The Arvesen et al. atlas is 
quite acceptable for our purposes. 

The second instrument to be used was the 
Large Cassegrain Spectrograph (LCS) at the 
2.7m telescope. This was a long slit CCD spec- 
trograph with various gratings. The bulk of the 
observations were obtained with a grating set 
to cover the spectral range from 3000 - 5700A 
at 7. 5 A resolution. Occasionally, we tilted the 
grating further to the red to observe other fea- 
tures or used a higher resolving power grating 
to study some bands in more detail. The slit 
was ~ 140arcsec long and could be oriented at 
any position angle by rotating the instrument on 
the back of the telescope. Each pixel on the de- 
tector covered a spatial interval of 1.28arcsec. 
The slit width was variable but was typically set 
to 2 arcsec for comet observations and 10 arcsec 
for standard star observations (so that no flux 
was lost). Generally, solar analog stars were 
observed at least once per run and often every 
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Figure 1: A spectrum of comet 122P/deVico obtained 
with the LCS on 27 September 1995 UT when the comet 
was at 0.69 AU from the Sun. The various molecular bands 
in our bandpass are marked on the spectrum. 

night. When not observed for whatever reason, 
convolution of the slit function with the Arvesen 
et al. (1969) atlas was used for a solar spectrum. 

With the LCS, we would always observe with 
the optocenter imaged on the chip at least once 
every night. We would then sometimes translate 
the slit under very accurate telescope control to 
extend our coverage into the coma. Some nights 
we also observed at more than one position an- 
gle. 

The third instrument was the Electronic Spec- 
trograph Number 2 (ES2), used on the 2.1m tele- 
scope. This instrument is very similar to the 
LCS and even used the same detectors. The slit 
position angle could also be changed by rotat- 
ing the instrument on the back of the telescope. 
Unlike the LCS, we could not image additional 
regions of the coma because the 2.1m telescope 
cannot be positioned accurately enough. 

Figure [1] shows a spectrum obtained with the 
LCS of comet 122P/de Vico on 27 September 
1995 UT. This is a very clean spectrum covering 
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1.28 X 2arcsec around the optocenter. The vari- 
ous features observed in our bandpass are noted 
on the figure. These are all molecular emission 
bands and are composed of many individual lines 
that we do not resolve. This comet had a very 
low dust-to-gas ratio so there is very little con- 
tinuum to be seen. The OH, NH and CN fea- 
tures tend to be relatively sharp and narrow and 
are easy to see and measure. The C3 and C2 
bands are broad and can easily be lost in noise 
when weaker, in spectra with lower signal/noise 
or when the continuum due to sunlight reflecting 
off the dust is very strong. The latter is a prob- 
lem because the continuum removal can repre- 
sent a large source of uncertainty and add noise 
to the result, even when the signal/noise of the 
comet spectrum is very high. The continuum de- 
creases faster with cometoccntric distance than 
does the C2 or C3 band strength, so this is more 
of a problem on the optocenter points than fur- 
ther into the coma. The CH is always very weak 
and difficult to measure. In the vast majority of 
the spectra reported in this paper, our bandpass 
included the features due to CN, C3, CH, and C2. 
We could never measure OH and NH2 simulta- 
neously because, as shown in the spectrum, the 
wavelength separation is beyond the bandpass of 
the instrument (only part of the NH2 is seen on 
this spectrum). In addition, high airmass makes 
observation of the OH very difficult (though iron- 
ically, the de Vico spectrum was obtained when 
the comet was at 2.5airmasses and we were able 
to observe the OH in this comet even when it 
was at 3.8 airmasses (!) because the comet was 
so close to the Sun). 

The IDS was used from 1980 through 1989. 
We started using the LCS in 1985. Thus, there 
was a period from 1985 - 1989 when we inter- 
changed which instrument was used (for reasons 
that have long since been forgotten!). The LCS 



was used with an RCA 512 x 512 pixel CCD 
from 1985 through 1986. This detector did not 

have UV response so was used in a redder set- 
ting but with more limited bandpass because of 
the smaller number of pixels than the subsequent 
CCD. From 1988, we used a Texas Instruments 
800 X 800 pixel CCD that was UV sensitized. We 
also used the Texas Instrument CCD on ES2. 
In all cases, the telescope was tracked at the 
comet's rate of motion on the sky and we guided 
the telescope directly off the slit. Thus, posi- 
tional accuracy was generally quite good. 

As would be expected with a survey that 
spanned so many years, we changed our observa- 
tional philosophy during the course of the survey. 
When we started, we wanted to be as inclusive 
as possible and would observe as many comets 
to V~19.5 as we could. We therefore sometimes 
observed comets at heliocentric distances greater 
than 2.5 AU, where they were generally not very 
active. We sometimes also sacrificed signal/noise 
for the opportunity to observe another comet. 
We might only observe a particular comet on a 
single night. 

As our survey progressed and as we changed 
instruments, we tended to concentrate more on 
targets that were brighter and at smaller helio- 
centric distances. Wc chose opportunities to ob- 
serve the same comet repeatedly over opportu- 
nities to observe additional comets. In addition, 
comets of special interest (mission targets, close 
approaches, etc.) began to dominate the sur- 
vey. Starting in 1995, whenever a comet was 
bright enough, we would observe at very high 
spectral resolution instead of with the instru- 
ments described here. Those observations are 
not included in this paper. 
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3. Reduction and Analysis 

The routine reduction of the data is very sim- 
ilar for each of the instruments. First, we re- 
moved the bias for the CCDs (the IDS did not 
have a bias) and flat fielded the data. Since 
the CCD instruments (LCS and ES2) are long 
slit instruments, we extracted the stellar spectra 
using variance weighting (an estimation of the 
variance created from a noise model based on 
the gain and rcadnoise parameters and a smooth 
profile function) of the slit profile to collapse the 
stellar spectral image into a ID spectrum. The 
comet spectra for the CCD instruments were left 
as 2D spectral images. For all instruments, we 
observed arc lamps in order to define the wave- 
length correspondence of each pixel. For the 
IDS, we used a fifth or sixth order polynomial to 
fit the dispersion. For the LCS, we needed to fit 
the wavelength in both the spectral and spatial 
directions. We used a fifth order polynomial in 
the spectral direction and a fourth order poly- 
nomial in the spatial direction. We would use 
our spatial fit to transform the spectral image so 
that the spectra were aligned with rows and the 
spectrum was no longer curved along the chip. 
Finally, the cometary and solar analog spectra 
were fiux calibrated with the standard star ob- 
servations to convert the spectra into flux versus 
wavelength. 

The cometary spectrum consists of an amal- 
gamation of two types of contributions. The 
coma gas yields a molecular band spectrum in 
emission, generally from resonance fluorescence. 
Sunlight reflecting off the dust creates a contin- 
uum underlying the gas spectrum. In addition, 
there is a third spectrum superimposed on the 
cometary spectrum and that is the spectrum of 
the Earth's atmosphere (the telluric night sky 
emission spectrum). 



For the purposes of this paper, we are inter- 
ested in the quantities of the gases that are im- 
aged through our aperture. Thus, we need to 
remove both the telluric night sky emission spec- 
trum and the solar spectrum. The most promi- 
nent telluric feature in our spectrum is the O (^S) 
feature at 5577A. The telluric spectrum also has 
bands in the UV that must be removed. To re- 
move the telluric spectrum, we would use either 
a dedicated sky observation (always the case for 
the IDS), obtained well removed from any comet, 
or we would use the spectrum at the end of the 
long slit (for the LCS and ES2) if a comet was 
not very extended and there were no gas fea- 
tures seen in that part of the spectral image. 
We would weight the sky spectrum by using the 
strength of the 5577A band above the local con- 
tinuum and then would subtract the sky spec- 
trum from the comet spectrum. Note that, for 
many comets, there is the C2 = — 1 band 
in the region of the 5577A line so we needed to 
define the continuum with respect to this band. 
Because the 5577A feature is by far the domi- 
nant telluric signal, slight mis-weighting of the 
telluric spectrum was not terribly critical except 
for the C2 Ai; = — 1 band. We never use this 
band in our analysis because there could be un- 
certainties induced in it by this process. We see 
no other strong telluric features (McDonald Ob- 
servatory does not suffer from street lighting in- 
duced features). There are some weak features 
seen, including the Herzberg bands of O2, pre- 
dominantly below 3600A. 

Some comets have very strong continuum 
spectra because their comae contain a great deal 
of dust to refiect sunlight; others have very weak 
continua. To remove the solar spectrum, it 
would be ideal to have observations of the Sun 
observed with the same instrument. This is not 
possible. The next best thing would be to ob- 
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serve a star that is a strong match to the prop- 
erties of the Sun. None of the tabulated solar 
analogs is fainter than V=10 so we were not able 
to observe them with the IDS, as noted above. 
For the LCS and ES2, we could observe solar 
analogs but considerations such as weather or 
availability meant that we did not always ob- 
serve one. The final resort for a solar spectrum 
was to use a full disk spectrum of the Sun. These 
are typically observed at high spectral resolution 
and must be convolved with a slit function to 
match our instrumental resolution. We adopted 
the solar spectrum of Arvesen et al. (1969) for 
this purpose. 

It should be noted that the small particles of 
dust in the coma result in the spectrum of the 
reflected sunlight being reddened when it is re- 
flected ofi' the dust. Thus, the solar spectrum 
does not exactly match the continuum we wish 
to remove. To fix this problem, we measured the 
flux in up to 11 continuum bandpasses in both 
the comet and the solar spectrum and weighted 
the solar spectrum to match the same color. Use 
of so many continuum regions allows us to match 
any non- linear reddening that exists. We used a 
smooth spline function to fit these continuum re- 
gions. The resultant weighting curve was never a 
straight line. This color weighted spectrum was 
then removed from the comet spectrum. This 
process would alleviate any issues with the UV- 
blue flux calibration of the Arvesen et al. (1969) 
atlas or with stars of slightly different colors than 
the Sun. 

In the case of the IDS, each observation would 
contain spectra of two positions - often the op- 
tocenter and 52 arcsec away. Each position 
was handled separately in the process described 
above. The off position still included the coma 
of the comet so was not used as a source of a 
sky spectrum unless the off^ spectrum was de- 



void of any emissions. For the LCS and ES2, 
the long slit covers >110 pixels spatially. For 
high signal/noise LCS or ES2 spectra, we would 
treat each spatial pixel separate spectrum 
for the telluric and solar removal. For lower sig- 
nal/noise observations, we would sometimes bin 
together pixels in the spatial direction and then 
apply the removal of the telluric and solar spec- 
tra (the concept of binning does not apply to the 
IDS). Figure [2] shows an example of handling 
the optocenter pixel for comet d'Arrest on 28 
September 1995. The optocenter spectrum has 
the strongest continuum. The continuum due 
to the dust weakens faster with spatial distance 
from the optocenter than does the gas emission. 

Once the telluric and solar spectra were re- 
moved, we could now measure the quantity of 
gas in each emission band. First, we fit a contin- 
uum to the region around each of the molecular 
bands of interest, removed the continuum and in- 
tegrated the fiux above the continuum. Though 
there should be no continuum left after removal 
of the telluric and solar spectra, in reality, the 
color weighting is imperfect because there are 
few true continuum regions for the weighting. 
The fitted continuum at this point provides a 
consistent level to integrate above. 

The integrated fluxes are then converted to 
column densities using the standard efficiency 
factors listed in Table [TJ These assume the coma 
is optically thin. This will always be the case in 
the outer coma. For a few of the most active 
comets, the collisional zone might be big enough 
to affect the inner pixel. However, our pixels are 
generally much larger than the affected region. 
One comet for which this would be a problem if it 
had been observed near 1 AU is comet C/1995 01 
(Hale-Bopp). However, as our observations were 
at 4au, this is not a problem for our data. 

Note that at this point we generally had col- 
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Table 1: Constants for Converting Fluxes to Column Densities 



Molecule 


Constant 
(ergs mol^"*^ sec^^) 


ReP 


OH 


f{Rh) 


1 


NH 


f{Rh) 


2 


CN 


f{Rh) 


3 


Cs 


12.42 


4 


CH 


12.98 


4 


C2 {Av = 1) 


12.62 


5 


C2 (Av = 0) 


12.35 


5 


NH2 


13.73 


6 



For OH, NH and CN 

log M = f{Rh) + log I. + 2 log Rh 
For C3, CH, C2, and NH2: 



log M = const + log L + 2 log Rh 
where M = column density, L = band intensity, R/i 
= heliocentric distance (au), Rh = heliocentric radial 
velocity (km/sec), f{Rh) is the fluorescence efficiency 
as a function of heliocentric distance. 
"References - 1: factor from Schleicher and 
A'Hearn 1982; 2: (/-factor from Kim et al. 1989; 3: 
Tatum and Gillespie (1977); 4: Cochran and Barker 
(1985); 5: Oliversen et al. (1985) and 17— factor from 
A'Hearn et al. (1985); 6: (7— factor from Tegler and 
Wyckoff (1989) 



umn densities for each observed molecule at 
multiple positions within the coma. For some 
comets, we might just have an optocenter point, 
while for others we had many hundreds of data 
points, covering distances to greater than 10^ km 
from the optocenter. Figure [3] shows an exam- 
ple of the column densities function of po- 
sition for comet Halley obtained using both the 
LCS and IDS. We needed a way to be able to in- 
tercompare comets that had disparate numbers 
of observed positions in the coma and a variety 
of observed molecules. We did this by model- 



ing the comets using the computationally simple 
Haser (1957) model. The Haser model is not 
very physical as it assumes that the cometary 
comae are spherically symmetric and that all 
gas only flows outwards from the nucleus. As 
we will show later, spherical symmetry is often 
violated. We also expect that dissociation will 
allow the daughter species to travel in any ran- 
dom direction. Despite these shortcomings, we 
still chose to use this model because it is an easy 
way to intercompare comets and to compare with 
other researchers. We adopted the Haser model 
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Table 2: Adopted Haser Model Scale Lengths 





Parent 


Daughter 






Scale length 


Scale length 




Molecule 


(km) 


(km) 


Reference" 


OH 


2.4x10* 


1.6x10^ 


1 


NH 


5.0x10* 


1.5x10^ 


2 


CN 


1.7x10* 


3.0x10^ 


3 


C3 


3.1x10^ 


1.5x10^ 


3 


CH 


7.8x10* 


4.8x10^ 


4 


02" 


2.5x10* 


1.2x10^ 


3 


NH2 


4.1 xlO^ 


6.2x10* 


5 



"^References - 1: Cochran and Schleicher (1993); 2: 
Randah et al. (1992); 3: Cochran (1986); 4: Cochran 
and Cochran (1990); 5: Cochran et al. (1992) 
C2 parent scales as R|'^, all others as R| 



scale lengths listed in Table [2] for this purpose. 
These are not necessarily the best fits to all of 
the data but offer a consistent set of parame- 
ters without too many degrees of freedom. Our 
choice of scale lengths, as compared with other 
authors, is discussed later. In addition, we use 
a modified Delsemme (1982) velocity scaling law 
of w = 0.85i?^'^ for the models. 

Table [3] is representative of the information 
that we gathered for each comet. For each comet 
and each night of observation, we list the obser- 
vational circumstances, the instrument, the po- 
sition angle of the slit for the LCS and ES2, mea- 
sured north through east (left blank for the IDS 
since the apertures were fixed on an east/west 
line) and the values of the measured produc- 
tion rates (denoted with a Q(X) where X is the 
species in question). For the production rates, 
we list the log of the production rate for each 
molecule along with the binning (for LCS and 
ES2) and the number of data points that went 
into the fit. Note that for some comets and 



some nights we would have more than one spec- 
tral image per night, resulting in many more 
data points than are covered by the slit in a sin- 
gle spectral image. These spectral images were 
sometimes, but not always, centered on different 
portions of the coma. 

The table here is for only a single comet 
(23P /Brorsen-Metcalf ) to give the reader an idea 
of the nature of the available data. As it is, this 
representation of the table must be sliced into 
three parts to fit the table onto a page for print- 
ing. Since there are over 650 lines of table with 
all of the comets, it is impractical to print all of 
the values. Instead, all of the data are available 
on-line in the supplementary materials for this 
paper. 

One of the columns in this table is labeled 
"Group" . This column lists the dynamical type 
of the comet: JFC = Jupiter Family comet (Tis- 
serand parameter > 2); HTC = Halley Type 
Comet (Tisserand parameter < 2 and P < 200 
years); LPC = Long Period Comet (P > 200 
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Table 3: Representative Production Rate Table (Broken into 3 parts) 



Comet 


Date 


PA 
(deg) 


Instr. 


(AU) 


R 

(km/sec) 


A 

(AU) 


Group 


Phot.? 


Log Q{OH) 
(mol / see) 


Bin 


Npts 


23P/Brorsen-Metcalf 


16 Jul 1989 


88 


LCS 


1.31 


-28.51 


0.84 


HTC 


n 


27.99 


2 


104 


09 Aug 1989 


175 


LCS 


0.90 


-29.68 


0.63 




y 










10 Aug 1989 


175 


LCS 


0.88 


-29.64 


0.63 




n 


27.98 


2 


54 




11 Aug 1989 


190 


LCS 


0.86 


-29.58 


0.63 




n 


28.60 


2 


225 



log Q(NH) 
(mol / see) 


Bin 


Npts 


log Q(CN) 
(mol/sec) 


Bin 


npts 


log Q(C3) 
(mol/see) 


Bin 


Npts 


Log Q(CH) 
(mol / see) 


Bin 


Npts 


25.37 


1 


148 


25.25 


1 


217 


24.64 


1 


184 


25.68 


1 


157 


26.11 


1 


467 


25.86 


1 


505 


25.11 


1 


449 


26.38 


1 


463 


26.13 


1 


175 




















26.17 


1 


503 


25.83 


1 


505 


25.19 


1 


389 


26.34 


1 


354 



Log Q(C2) 
(mol/sec) 
Av = 1 


Bin 


Npts 


Log Q(C2) 
(mol / see) 
Au = 


Bin 


Npts 


Log Q(NH2) 

(mol/scc) 


Bin Npts 


25.47 


1 


213 


25.46 


1 


217 






26.10 


1 


490 


26.09 


1 


496 


25.35 


1 277 


26.10 


1 


498 


26.11 


1 


504 







years). The next column is labeled "Phot?". 
The values for this column are either "y" or "n" 
to indicate whether or not it was photometric. 
If it was not photometric, then the production 
rates are not absolute. However, since all the 
wavelengths were observed simultaneously, the 
relative productions of species are still accurate 
since most clouds are gray (Wing 1967). For 
the LCS and ES2, we would sometimes obtain 
multiple spectral images even if not photomet- 
ric. For purposes of the analysis, we would scale 
all the column densities (as a group from a sin- 
gle spectral image) to the spectral image that 
was the brightest. For this scaling, all molecules 
were scaled by the same amount. That we could 
use a single scaling factor to bring molecules 
at disparate wavelengths into good alignment is 
another indicator that the clouds are generally 
gray. The resultant production rates are still 
not absolute but, again, the relative production 
of species is correct. 

When we observed at more than one slit po- 
sition angle with the LCS or ES2, we would 
keep the different position angles separate for the 
analysis. The position angle should not affect 



the production rates (as long as the optocenter 
is in each image) so this serves as a check on 
the derived production rates. The exception to 
this policy was the Encke observations of Octo- 
ber 2003 as we used so many different position 
angles that it was impractical to keep them sep- 
arate. 

Inspection of Table [3] shows that some of the 
entries for production rates are blank. In the 
case of Brorsen-Metcalf, these blanks are be- 
cause the bandpass did not include some of the 
molecules on some of the nights. On 16 July 1989 
and 11 August 1989, we observed with our nom- 
inal setup that covered OH but not NH2. On, 
9 August 1989, we changed the grating angle to 
observe NH2 but not OH. On 10 August 1989, we 
employed a higher resolving power grating and 
only observed OH and NH. 

For other comets, there are blank entries even 
when we covered a molecule in our setup. This 
means that we did not detect the molecule, gen- 
erally because the data were too noisy. We con- 
sider that we should be able to measure CN, 
C3, and C2 in all comets (OH and NH are af- 
fected by airmass and low detector response, CH 
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Flux calibrated optocenter spectrum 
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Sky spectrum removed 
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Solar continuum removed 
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Wavelength (A) 

Figure 2: This figure demonstrates the process of removal 
of the telluric night sky emission and solar spectra from 
the original spectrum of comet 6P/d' Arrest. The top 
panel shows the flux and wavelength calibrated spectrum 
of the comet before removal of either of these spectra. 
Note the strong feature at 5577A. This is the O (^S) band 
in the telluric spectrum. The O2 Herzberg bands are also 
marked. The continuum results from sunlight reflecting 
off the dust. The middle panel shows the same spectrum 
but with the telluric spectrum removed using an inde- 
pendent observation of the sky. The bottom panel shows 
that spectrum after the removal of the solar continuum. 
In this case, the spectrum of solar analog Hyades 64 was 
removed after weighting the continuum to match the col- 
ors of the star and the dust (see text). 



is very weak and NH2 is often off the chip). 
We can put upper limits on the column den- 
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Figure 3: The column densities as a function of position 
in the coma are shown for comet Halley. In this figure, we 
show the data sets from both the IDS (13 March 1986) 
and LCS (24 March 1986). For the IDS, we only observed 
4 bands, while we observed 6 bands with the LCS. Note 
the much larger number of data points obtained with the 
LCS, the result of the long slit. Superimposed on the 
data are the Haser model fits to the data. The production 
rates, log Q, are in mol sec~^; the velocity, vel, is in km 
sec~^; the heliocentric distance, Rh, in in AU. 



sities of these three molecules by looking at 
how the data in the bandpass of the emission 
band deviates from the continuum. To mea- 
sure an upper limit, we first fit a continuum 
to the first and last three points in the band- 
pass. Next, we computed a standard deviation, 
a = sqrt[T,{observed — continuum^ / {N — 1)], 
where N is the number of data points in the 
bandpass. The upper limit is then defined as 
U.L. = 0.75aB, where B = the bandpass in A. 
This is a 3a upper limit. 

Figure S] shows an example of a comet for 
which we detected the CN band but there is 
no obvious C3 or C2- This spectrum, of comet 
192P/Shoemaker-Levy at of 1.84 AU, has had 
the telluric and solar spectra removed. The 
bandpasses for the C3 and C2 Av = bands are 
marked on the figure. Included are the measured 
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Figure 4: A spectrum of comet 192P/Shoemaker-Levy is 
shown. CN is obviously visible but there is no sign of 
the C3 or C2. By measuring the noise in the bandpass 
where we would expect to detect these features, we can 
derive upper limits for the amount of gas that could be 
attributable to these species. These upper limits are 3a" 
limits as discussed in the text. 



integrated band intensity for CN along with the 
3(7 upper hmit integrated fluxes for C3 and C2. 
In this example, it should be obvious that there 
cannot be much C2 or C3 hidden in this spec- 
trum. 

In general, the upper limits are not very con- 
straining. If we observed all three of CN, C3, 
and C2 on at least one of the nights we observed 
a particular comet, we did not tabulate the up- 
per limits for other nights because we believed 
that the actual detections were more meaningful 
(the upper limit measurements are much more 
labor intensive than the usual column density 
measurements because we never developed the 
same degree of automation in their computation 
that we have for column densities). However, we 
have included upper limits for the cases when CN 
was detected but C2 and C3 were not both de- 
tected. We have only one case, comet C/1984K1 
(Shoemaker), when we detected CN and C3 but 



not C2 and one case, C/1992 Wl (Oshita), where 
we observed CN and C2 but not C3. 

Of our sample of 130 comets, there are 20 
comets for which we observed no molecules, 
not even CN. These are tabulated in Table HJ 
where we list the date, observing circumstances, 
group and the upper limits for the CN pro- 
duction rates. We do not include any other 
upper limits since ratios of upper limits are 
not at all meaningful. Inspection of Table H] 
shows that many of these comets were observed 
at large heliocentric distance (> 2.5 Au), where 
they would not be expected to produce much 
emission. However, there are a few exam- 
ples (108P/Ciffreo, 96P/Machholz 1, 7P/Pons- 
Winnecke, and 125P/Spacewatch) with observa- 
tions at < 1.75 AU, where we would have ex- 
pected to detect at least CN. These comets must 
be very low producers, indeed. 

Before we can compare comets, it is impor- 
tant to understand the uncertainties in our mea- 
surements. As with all data, these include both 
systematic and random errors. One source of 
the systematic errors is our knowledge of the g- 
factors. We differ from A'Hearn et al. (1995) for 
the C3 and CN ^(-factors. The different (7-factors 
will make it hard to intercompare with other ob- 
servers but should not affect intercomparisons of 
comets within a single data set. 

For C3, no good oscillator strengths have been 
measured in the laboratory since C3 is not a sta- 
ble molecule. Jorgensen et al. (1989) reviewed 
the laboratory measurements and found a range 
in measured oscillator strengths of a factor of 
6. We long ago adopted a value of / = 0.016 
(Clegg and Lambert 1982) while A'Hearn et al. 
use 0.001. Roueff et al.. (2002) adopted a value 
of / = 0.0146 whereas Adamkovics et al. (2003) 
adopted / = 0.016, in agreement with our value. 
Our value is consistent with the current liter- 
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Table 4: Comets With No Molecules Detected 



Comet Date Instr. R,, R A Group Log Q(CN) 

(AU) (km/sec) (AU) Upper Limit 

(mol/sec) 



C/1980 El (Bowell) 



C/1981 HI (Bus) 



16 Oct 


1985 


IDS 


2, 


.41 


-3.52 


2, 


.16 


JFC 


<24, 


.04 


16 Nov 


1985 


IDS 


2, 


.36 


-2.54 


2, 


.47 




<24, 


.07 




1 Q80 


IDS 


o ■ 


38 


-1 1 8^ 

XX .(JO 


U . 


^6 


T,pr! 

XJX 


^9^ 






1 Q80 


IDS 

XX-/ kJ 


o ■ 


36 


-1 1 83 

XX. oo 




4Q 






9Q 


Tan 


1 Q81 


IDS 

XX-/ o 




93 


-1 1 6^ 

XX .KJO 




07 






1 Q 

. X t/ 


07 Feb 




IDS 

xxy kJ 


5 


on 

• \J\J 


-11.29 


4 


33 




<25 


71 

• 1 X 


1 1 Fph 


1 981 


IDS 


4, 


Q8 


-11.25 


4, 


96 




^25 


56 


07 Jun 


1981 


IDS 


4, 


.27 


-9.72 


3, 


.88 




<24, 


.91 


24 Jan 


1982 


IDS 


3, 


.39 


-2.24 


3, 


.92 




<25, 


.67 


25 Jan 


1982 


IDS 


3, 


.39 


-2.19 


3 


.90 




<25, 


.41 


20 Apr 


1982 


IDS 


3, 


,38 


1.84 


2, 


.82 




<24, 


.90 


18 May 


1982 


IDS 


3, 


,42 


3.12 


2, 


.57 




<25, 


.25 


26 Jun 


1982 


IDS 


3. 


.51 


4.77 


2, 


.50 




<25, 


.40 


16 Jul 


1982 


IDS 


3, 


.57 


5.55 


2, 


.61 




<25, 


.21 


17 Jul 


1982 


IDS 


3, 


,58 


5.59 


2, 


.62 




<25, 


.17 


18 Jul 


1982 


IDS 


3, 


,58 


5.63 


2, 


.63 




<25, 


.11 


17 Aug 


1982 


IDS 


3. 


.69 


6.68 


2, 


.99 




<25, 


.14 


15 Sep 


1982 


IDS 


3, 


.81 


7.58 


3 


.47 




<25, 


.15 


13 Oct 


1982 


IDS 


3, 


,94 


8.33 


3, 


.99 




<25, 


.06 


10 Jun 


1983 


IDS 


5, 


,36 


11.50 


5, 


.01 




<25, 


.33 


13 Jun 


1983 


IDS 


5, 


.38 


11.52 


4, 


.99 




<25, 


.08 


07 Jul 


1983 


IDS 


5. 


.55 


11.63 


4, 


.83 




<25, 


.30 


05 Aug 


1983 


IDS 


5. 


.74 


11.73 


4, 


.78 




<25, 


.07 


10 Sep 


1983 


IDS 


5, 


.99 


11.83 


5, 


.02 




<25, 


.10 


05 Oct 


1983 


IDS 


6, 


,16 


11.88 


5, 


.41 




<25, 


.11 


06 Jun 


1981 


IDS 


2, 


,53 


-4.46 


1, 


.91 


LPC 


<23, 


.73 


07 Jun 


1981 


IDS 


2. 


,52 


-4.39 


1, 


.92 




<23, 


.63 
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Table H Comets With No Molecules Detected (cont.) 



Comet Date Instr. Rh R A Group Log Q(CN) 









(AU) (km/sec) 


(AU) 




Upper Limit 

I iilVjl / OC^j J 


108P/Ciffreo 


12 Dec 1985 


IDS 


1.75 


3.94 


0.79 


JFC 


<23.32 




ni Till 1 QSQ 




1.98 


-7.46 


1 1 Q 






C/1983 Nl (IRAS) 


02 Apr 1984 


IDS 






3.27 


LPC 


<24.49 


48P/ Johnson 


08 Jul 1983 

06 Aug 1983 

07 AiiP 1 Q83 


IDS 
IDS 
IDS 


2.44 
2.44 


—4:. "4:0 

-3.78 
-3.76 


1.50 
1.58 
1 59 

± .OiJ 


JFC 


<23.61 
<23.80 
^93 91 


C/1999 S4 (LINEAR) 


25 Jan 2000 

26 Jan 2000 


LCS 
LCS 


2.97 
2.96 


-21.07 
-21.10 


2.80 
2.81 


LPC 


<24.97 
<24.45 


77P /Longmore 


04 Jan 1981 
17 Apr 1988 


IDS 
IDS 


2.98 
2.65 


-5.53 
-4.39 


2.31 
2.09 


JFC 


<23.17 
<23.23 


96P/Machholz 1 


13Jun 2007 
14Jun 2007 


LCS 
LCS 


1.62 
1.64 


26.78 
26.62 


0.67 
0.68 


JFC 


<23.68 
<23.43 


7P /Pons-Winnecke 


06 Apr 1989 
09 May 1989 


IDS 
IDS 


1.98 
1.73 


-13.12 
-12.56 


1.19 
1.17 


JFC 


<23.55 
<23.76 


30P/Reinmuth 1 


18 Nov 1987 

19 Nov 1987 
24 Dec 1987 


IDS 
IDS 
IDS 


2.38 
2.38 
2.22 


-8.35 
-8.34 
-7.57 


1.43 
1.42 
1.31 


JFC 


<24.26 
<23.81 
<24.42 


94P/Russell 4 


04Jun 1984 


IDS 


2.38 


5.24 


1.70 


JFC 


<24.25 


C/1983 Rl (Shoemaker) 


05 Oct 1983 

06 Dec 1983 


IDS 
IDS 


3.38 
3.35 


-2.25 
0.56 


2.45 
3.46 


LPC 


<24.39 
<24.48 
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Table H Comets With No Molecules Detected (cont.) 



Comet Date 


illbLl . 




R 


A 
LA 


Group Log Q(CN) 








(AU) (km/sec) 


( ATT'i 




Upper Limit 
















(mol/sec) 


C/1989 A5 (Shoemaker) 07 Feb 1989 




2.65 


-1.39 


i. ( u 


LPC 


<25.17 


08 Feb 


1989 


IDS 


2.65 


-1.32 


1.70 




<24.70 


04 Apr 1989 


IDS 


2.67 


2.66 


2.43 




<24.04 


08 May 1989 


IDS 


2.74 


4.93 


3.00 




<24.68 


128P/Shoemaker-Holt 1 18 Nov 1987 


TFtQ 


o.ZZ 


-z.yb 


9 /in 


Jr U 


<z4.oo 


11 Dec 


1988 


TDQ 
iUo 


3.25 


3.15 






^9/1 /I Q 


121P/Shoemaker-Holt 2 05 Apr 1989 


IDS 


3.01 


4.56 


2.27 


JFC 


<24.15 


56P/Slaughter-Burnham 25 Aug 


1981 


IDS 


2.63 


-3.25 


1.87 


J r 


^90 on 


26 Aug 


1981 


lUkj 


2.63 


-3.22 


1 .ou 




<23.36 


28 Aug 


1981 


IDS 


2.62 


-3.15 


1.84 




<23.56 


27 Sep 


1981 


IDS 


2.58 


-2.07 


1.61 




<23.50 


28 Sep 


1981 


IDS 


2.58 


-2.04 


1.60 




<23.47 


24 Oct 


1981 


IDS 


2.55 


-1.03 


1.58 




<23.81 


23 Jan 


1982 


IDS 


2.59 


2.54 


2.51 




<23.81 


25 Jan 


1982 


IDS 


2.60 


2.58 


2.52 




<23.45 


74P/Smirnova-Chernykh 26 Oct 


1981 


IDS 


4.34 


-1.98 


3.44 


JFC 


<24.27 


13 Oct 


1982 


IDS 


3.91 


-2.07 


3.79 




<24.32 


14 Oct 


1982 


IDS 


3.91 


-2.07 


3.78 




<23.98 


16 Dec 


1982 


IDS 


3.84 


-1.95 


2.92 




<24.03 


06 Dec 


1983 


IDS 


3.57 


-0.45 


3.58 




<24.21 


08 Dec 


1983 


IDS 


3.57 


-0.44 


3.55 




<24.14 


25 Feb 


1984 


IDS 


3.56 


0.02 


2.61 




<23.93 


31 Mar 1984 


IDS 


3.56 


0.23 


2.62 




<24.46 


125P/Spacewatch 15Jun 


1996 


LCS 


1.57 


-3.12 


1.08 


JFC 


<23.65 
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ature. However, the oscillator strength is, in 
essence, a scaling factor for the column densities, 
so our C3 column densities would be different by 
a constant factor from those computed with a 
different oscillator strength. 

CN is a bit more complicated than C3 because 
one needs to account for the Swings effect. Thus, 
the fluorescence efficiency is a function of the he- 
liocentric radial velocity of the comet. We have 
adopted the older Tatum and Gillespie (1977) 
calculation of the fluorescence efficiency while 
A'Hearn et al. use a formulation that factors 
in the heliocentric distance of the comet. Recent 
calculations by Schleicher (2010) are more com- 
plete, covering more heliocentric distances than 
previous calculations. Schleicher concludes that 
the values at a particular heliocentric distance 
vary by about a factor of two with velocity but 
that the trend of the fluorescence efficiency with 
radial velocity at different heliocentric distances 
can be radically different. Schleicher's calcula- 
tions were published after we had completed the 
reduction and analysis of much of our data. We 
did not switch fluorescence efficiencies because of 
the large investment in manpower at that point. 
The effect of using the older calculation is not 
a simple constant offset, as with the C3 oscilla- 
tor strength. It will be less than the full range 
of variation for the fluorescence efficiency, or a 
factor of 2 (0.3 in the log), and will be responsi- 
ble for increasing the error bars when values are 
averaged later in this paper. 

Another source of systematic errors is our 
knowledge of the scale lengths for the model. 
We do not take into account any change in 
the strength of the solar UV spectrum nor do 
we take into account any multi-step dissocia- 
tion processes. Inspection of the model fits in 
Figure [3] show examples of good fits (CN and 
C3) and poor fits (C2 and NH2). The C2 is al- 



most always a bad fit in the very inner coma 
because C2 is probably a granddaughter prod- 
uct of its source and the Haser model does not 
account for the two-step production. In addi- 
tion, some of these scale lengths work for some, 
but not all, comets. This problem was noted by 
LSll who chose to fit the scale lengths for each 
comet. We do not do this because then there 
are too many free parameters in the intercom- 
parison of comets. Indeed, if the UV insolation 
does not change, the scale lengths are, to first 
approximation, determined entirely by the solar 
flux and the distance of the object from the Sun. 
While collisions might affect the inner coma, the 
vast majority of the coma is unaffected. Our 
flts rarely end up being dependent on the in- 
ner coma data points, anyway, since we have so 
many data points at larger cometocentric dis- 
tances. We described the process of defining the 
scale lengths in Cochran (1986) and use these or 
updated scale lengths uniformly for all of our ob- 
servations. These scale lengths were defined us- 
ing observations from five comets observed from 
0.74 to 1.81 AU. We include an expanded discus- 
sion of the scale lengths later when we intercom- 
pare different data sets. 

Probably the most important source of sys- 
tematic error is the Haser model itself. This 
model assumes that the cometary coma is spher- 
ically symmetric and that gas can only fiow out- 
wards. With our long slit observations, we ob- 
serve that spherical symmetry is rarely a good 
approximation. The effect of using the Haser 
model is, in affect, to azimuthally average the 
coma observations and assume this is a reason- 
able proxy for the production. Figure [5] is an 
(albeit extreme) example of the types of asym- 
metry we observe. In these observations of comet 
2P/Encke, obtained on 24 October 2003, we ro- 
tated the slit to 7 different position angles. Since 
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Figure 5: The CN column densities for 2P/Encke ob- 
served on 24 October 2003 arc shown. Each panel rep- 
resents a different slit position angle (and that position 
angle ±180° on the other side of the optocenter). There 
are two spectral images at each position angle. Clearly, 
comet Encke's coma is not symmetric. 



the comet was somewhere near the center of the 
detector, spatially, we actually obtained data on 
two sides of the optocenter (the listed position 
angle and the position angle ± 180°) for each 
listed position angle. Additionally, at each posi- 
tion angle we obtained two spectral images (that 
are shown with different symbols in each panel) . 
Inspection of this figure shows that Encke's coma 
is far from symmetric! The asymmetries ob- 
served are not due to clouds or other errors. This 



is demonstrated by the fact that both spectral 
images at a single position angle are in excel- 
lent agreement (no scaling was done). At some 
position angles (e.g. 70°, when the slit was po- 
sitioned perpendicular to the Sun/tail line), the 
coma looks reasonably symmetric, while at oth- 
ers (e.g. 160°, along the Sun/tail line), there is 
a factor of 10 difference in column density by 
the time we are a few thousand km from the op- 
tocenter. A simple model like the Haser model 
cannot do these data justice. A more compli- 
cated model, such as proposed by Ihalawela et 
al. (2011), allows for deviations from spherical 
symmetry. However, it should be noted that if 
we independently fit the different position an- 
gles in these Encke data with the Haser model 
(using both sides of the optocenter together in 
a single fit), we get the same CN production 
rate for each position angle (even though the fits 
look bad). Asymmetries are observed in many 
comets. Combi and Fink (1993) showed that 
Haser model scale lengths could be derived to ac- 
commodate time variable asymmetries in comet 
Halley. 

Another source of potential systematic error 
is the flux calibration. The fluxes of the stan- 
dard stars are probably only good to 5 10%. In 
addition, though we used a large aperture for ob- 
servations of the standard stars, there could be 
a loss of light from the aperture on some nights. 
This would cause a decrease in the amount of 
light we observed for the star and our sensitivity 
function would then cause a slightly higher flux 
for the cometary spectra. This would be an ex- 
tremely rare occurrence, especially for the LCS, 
since we used a lOarcsec wide slit. 

Clouds would also cause a diminution of the 
flux. They can affect the comet and standard 
stars by differing amounts. Thus, it could cause 
the cometary flux to appear brighter or dimmer. 
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As noted above, clouds are gray and we observe 
all wavelengths at the same time, so we can in- 
tercompare the different bands in a ratio sense. 
In addition, with the long slit, we observe many 
different positions at the same time, so the col- 
umn densities as a function of position maintain 
an accurate relationship. It is only when we were 
obtaining more than one spectral image that we 
had to scale the data. As noted above, we would 
use the same scale factor at all wavelengths. 

An additional source of error was continuum 
removal. This is required in order to study 
the gas. Sometimes the noise in the continuum 
meant that the process was not as accurate as 
we would desire. We mitigate this uncertainty 
in large part by fitting a continuum under the 
band prior to computing the integrated band in- 
tensity. This corrects for any continuum removal 
errors. 

Finally, an important stochastic error is the 
photon statistics. For some of our observations, 
such as the Encke data shown in Figure [5l it is 
clear that the signal/noise was very high. For 
others, the signal/noise is not at all good. This 
is especially true for some of the IDS data where 
we chose to obtain spectra of another object in- 
stead of building signal on the object we were 
observing. As is normal for stochastic errors in 
data, we can lower the impact of this noise by 
collecting more data, i.e. by observing more po- 
sitions in the coma. For our long slit CCD data, 
this is a natural consequence of the long slit. 

We can judge the affect of the noise on our 
answer by using a Monte Carlo model to simu- 
late changing the noise. Photon statistical errors 
affect each point individually. Thus, for a given 
signal, the measured value could vary based on 
the signal/noise. To simulate the effect, we took 
a data set and we changed the value of each data 
point by selecting randomly from a set of normal 
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Figure 6: The effects on the column densities due to pho- 
ton noise are illustrated in these panels. See the text for a 
discussion of the model. The top panels are three different 
C2 data sets with 433 individual column densities (Halley 
24 March), 30 column densities (Stephan-Oterma 9 De- 
cember) and only 7 column densities (Halley 13 March). 
In the top panels, the original data are shown with the 
best fit Haser model. In the bottom panels, the data are 
shown after we draw new values from a normalized Gaus- 
sian distribution for each point based on 30% potential 
error. The model fit is the same one as the real data 
and looks good for all three data sets. However, while 
the same model is correct for the Halley 24 March and 
Stephan-Oterma data, the Halley 13 March data actually 
required lowering the log C2 column density by 0.05. This 
change in best fit is based on 1000 simulations of noise. 



Gaussian deviates. Once each data point was 
modified, we fit the Haser model and recorded 
the value for the production rate. We then did 
this same process repeatedly. Finally, we figured 
the mean and standard deviations of the derived 
column densities (in linear space) with the model 
that was fit to the real data. We tested the er- 
rors by running the simulation with 1000 sets of 
modified data. The range of values that can be 
selected are generally set by the signal/noise of 
the data. Figure [6] shows examples of the pro- 
cess in three different C2 data sets - one with 
433 data points (Halley from 24 March 1986), 
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one with 30 data points (Stephan-Oterma from 
9 December 1980) and one with 7 data points 
(Halley from 13 March 1986). In ah three cases, 
we added 30% noise, though it is apparent from 
the scatter in the data that 30% noise is an over 
estimate, at least for the 24 March Hahey (ac- 
tual noise about 2.5%) and for Stephan-Oterma 
(noise of about 7%). As shown in the figure, even 
with the added noise, the Haser model gives the 
same values for the production rates for those 
two data sets. For the 13 March Halley data, 
30% error would result in a change of the pro- 
duction rate of only 0.05 in the log. Thus, to 
the extent that the model and the scale lengths 
are appropriate, our column densities are very 
robust when we have more than 25 data points. 
Fewer data points means that the uncertainty in 
the production rates might be 5 - 20% (stochas- 
tically). This is why the number of data points 
is included in the table that is available on line 
and shown as an example in Table [3l 

In addition to looking at the number of data 
points and how much scatter appears in plots of 
column density versus position, one can check 
consistency using C2 observations. We measure 
both the C2 Af = and Av = 1 bands in our 
spectra. The Av = band has a larger inte- 
grated band intensity than the Av = 1 band. 
This is accounted for by the different constants 
used for conversion of the intensities to the col- 
umn densities. Therefore, the two bands should 
give independent measurements of the C2 column 
densities. Figure [7] shows our derived column 
densities for all comets on all nights on which 
we measured both bands. There should be a 1:1 
correlation. Inspection of Figure [7] shows that 
in general this is true, though there are some 
outliers. We use different symbols for the IDS 
and LCS data. Most of the outliers are the 
IDS data, as would be expected from the much 
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Figure 7: The measured values of the production rates 
for the C2 An = 1 and C2 Av — bands are shown plot- 
ted against one another. We would expect to derive the 
same production rates from either band. This is what we 
observe. The IDS data are denoted with squares while 
the LCS data are denoted with triangles. The uncer- 
tainties in the IDS production rates are generally higher 
than for the LCS production rates since we observe so 
many more positions with the LCS. The disagreement of 
the production rates increases at the smallest production 
rates, indicative of larger uncertainties for weaker comets. 



smaller numbers of points obtained for the IDS, 
resulting in larger error bars on the derived pro- 
duction rates. In particular, we see the greatest 
deviations from the correlation at the smallest 
values for the column densities. This is not sur- 
prising as the weakest features are the ones with 
the largest uncertainties. When log Q for either 
band is >25, the agreement is generally quite 
good. A few outliers with log Q > 25 exist, and 
were checked, but there was nothing to pick one 
value over the other. 



18 








26 27 26 23 24 25 26 27 

log Q(CN) [mol sec"'] 



23 24 25 26 27 



Figure 8: This figure illustrates the observed production 
rate trends. In each panel, we show the log of the produc- 
tion rate of a particular molecule plotted against the log of 
the CN production rate. A slope equal 1 line is included 
on the plot to guide the eye (these lines are not fits). 
We see that the production rates of all of the molecules 
correlate well with the CN production rates. Except that 
JFCs are generally lower gas producers, we see no obvious 
trends with dynamical type. 



4. Trends in the data 

Figure El shows the observed column densities 
for ah comets on each night they were observed, 
with OH, NH, C3, CH, C2 and NH2 being plotted 
against CN. While OH, a daughter of the most 
abundant ice H2O, would be the most desirable 
molecule to use to compare one comet against 
another, we chose to plot abundances versus CN 
because CN is the most commonly observed fea- 
ture in our spectra. Inspection of this plot shows 
that the production rates of each species, includ- 
ing OH, are well correlated with the production 
rate of CN. Thus, comparing to CN is a rea- 
sonable approach. In this figure, we denote the 
different dynamical types of comets by different 
symbols. We see no obvious trend with dynam- 
ical type except that JFCs tend to produce less 
of all species. HTCs and LPCs seem to produce 



comparable amounts of gas. 

Eighteen of the comets were observed on only 
one night each. An additional 20 comets were 
only observed on two nights each. The remaining 
72 comets were each observed on at least three 
different nights or position angles, with comet 
9P/Tempel 1 having 45 observations. Clearly 
then, some comets are dominating the trends in 
Figure El We can mitigate this effect by combin- 
ing all the data for an individual comet. 

In order to combine the data, we first took 
ratios of the production rates of each species 
with respect to CN, for each date when both 
molecules were detected. This had the effect 
of removing any weather-related differences be- 
tween comets, as clouds would afi^ect all species 
equally. In addition, unless a comet is inhomo- 
geneous, we expect that all the species should 
be produced in the same ratio at different helio- 
centric distances. Figure shows a comparison 
of the change in production rates function 
of heliocentric distance versus the change in pro- 
duction rate ratios as a function of heliocentric 
distance for comet IP/Halley. The production 
rates show a dramatic increase when the comet 
is at smaller heliocentric distances. In contrast, 
the production rate ratios are constant with he- 
liocentric distance. Fink (1994) also found a con- 
stant mixing ratio with heliocentric distance for 
comet Halley. 

Figure [TU] shows the observed production rate 
ratios as a function of log heliocentric distance 
(using the log for the heliocentric distance just 
spreads out the data along the x axis since the 
vast majority of data were obtained between 1 
and 2.5 Au). With the possible exception of C3, 
there are not any obvious trends with heliocen- 
tric distance, as would be expected if the species 
are outgassing together. We also see no obvious 
trends with dynamical type. The dashed lines in 
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Figure 9: This figure shows how the gas production 
rates change as a function of heliocentric distance (up- 
per panels) for comet Halley. We see each species show- 
ing increased production at smaller heliocentric distances. 
Closed symbols denote photometric nights while open 
symbols mean it was not photometric. In the lower panels 
we plot the ratio of the production rate of each species to 
that of CN, with the y axis having the same overall range 
(3.8 dex) as in the upper panels. Notice that the produc- 
tion rate ratios are constant with heliocentric distance. 

the C2 and C3 panels denote our delimiters for 
depleted comets, discussed below. 

The C3 ratios for log heliocentric distance 
> 0.5 appear to be systematically higher than 
the rest. There are 8 data points at these large 
distances but the data represent only 5 comets 
(there are three Hale-Bopp points). One of those 
is comet IP/Halley. Inspection of Figure [9] sug- 
gests this point is consistent with other Halley 
observations at smaller heliocentric distances. 
Comet C/1983 01 (Cernis) has the highest C3 
ratio at this distance. This is the only C3 value 
measured for Cernis so we do not know if the 
comet is different or the value is high because of 
the large heliocentric distance. Comet C/1993 
Al (Mueller) has the next highest value but is 
also on the high side at 2 AU {log distance = 
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Figure 10: The production rate ratios are plotted as a 
function of log heliocentric distance for all instances where 
we observed, ratioed to CN for the same night. Each point 
represents one comet on one particular night. The differ- 
ent dynamical types are denoted with different symbols 
with the symbols defined in the OH panel. There are 
no obvious trends with heliocentric distance. The dashed 
lines for C2 and C3 represent the dividing line between 
carbon-chain typical and depleted comets (see text). 

0.31). Comet C/1996 PI (Wilson) has a high 
value of C3 at 3.5 AU but normal-to-low values 
at 1.34 and 3.13 AU, suggesting the high value 
at 3.5 AU may be in error. Three of the points 
belong to C/1995 01 (Hale-Bopp). We only ob- 
served this comet at large heliocentric distances 
as part of this project so we do not know if the 
comet is different or there is a heliocentric dis- 
tance effect. 

The constancy of production rate ratios with 
heliocentric distance means that we can compute 
average values for the ratios that are representa- 
tive of each comet at all observed heliocentric 
distances. The averages make the data for a 
comet with 3 observations as meaningful as a 
comet with 40 observations. In addition, since all 
species were observed simultaneously, the effects 
of clouds are canceled. Similarly, we have com- 
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Figure 11: Twelve comets were observed at a large enough 
range of heliocentric distances that we could investigate 
how the CN production rates scale with heliocentric dis- 
tance. We did not attempt to differentiate between pho- 
tometric and non-photometric conditions, nor did we sep- 
arate apparitions. These panels show the data along with 
the line that best fits the data (solid line). Also shown is 
a 1/R^ line (dotted line). The derived average CN val- 
ues for I'AU using the fit and the l/R^, law are shown in 
Table m 



dotted line shows a 1/R^ line. For some comets, 
there is a big difference between the fit and the 
1/R^ line. Table [5] lists the derived slopes and 
the averages from using a 1/R| scaling versus 
the fitted scaling. Not surprisingly, we get dif- 
ferent 1 AU CN values depending on what we as- 
sume. However, as we only use the 1 AU CN for 
Figure [T2] (described below) and that value does 
not affect our conclusions, the exact scaling used 
is irrelevant, as it only moves points left or right 
and we care most about vertical changes. None 
of the other species were scaled to 1 AU; all ratios 
were taken without scaling species to 1 AU. 
These average values are given in Table [6l 



puted an average CN production for each comet 
by computing the production the comet would 
have at 1 AU assuming a 1 /R| dependence of the 
production. 

For the discussion that follows, it is not critical 
if the CN varies differently than we have assumed 
because any deviation from this law will result 
in more scatter in the 1 AU values, and there- 
fore larger error bars. The exact values are not 
critical to our trends. For most of our comets 
we did not have a large enough range of helio- 
centric distance to investigate the correct scaling 
with heliocentric distance; for twelve comets, we 
do cover large heliocentric distance ranges and 
these are shown in Figure [TTl We fit a line to 
the data but did not in any way try to correct 
for weather or for changes from one apparition 
to the next. The solid line shows the fit while the 
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Table 5: The Heliocentric Dependence of the CN Production Rates 



Heliocentric Avg. Q(CN) Avg. Q(CN) 

Distance scaling at 1 AU at 1 Au 

Comet Range (au) Npts slope (fitted) (l/R-^) 
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Table 6: Average Production Rate Ratios 



log Q(CN) 



Average log production rate ratios 
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0.19 ± 0.164 


14 


0.06 ± 0.103 


20 


-0.11 ± 0.178 


9 


144P/Kushida 


26.27 ± 0.16 


4 


2, 


.01 ± 0.177 


4 


0.28 ± 0.014 


2 


-0.61 ± 0.077 


4 


0.20 ± 0.122 


2 


0.13 ± 0.038 


4 






C/1987 Al (Levy) 


24.74 ± 0.39 


3 












-0.61 ± 0.018 


3 


0.01 ± 0.190 


2 


0.14 ± 0.071 


3 


-0.02 ± 0.057 


2 


C/1990 Kl (Levy) 


26.36 ± 0.20 


5 


2, 


,76 ± 0.066 


3 


0.60 ± 0.098 


4 


-0.46 ± 0.142 


5 


0.32 ± 0.089 


4 


0.32 ± 0.075 


5 






C/1984 VI (LR) 


25.21 ± 0.24 


6 












-0.69 ± 0.087 


6 


0.13 ± 0.287 


5 


0.06 ± 0.056 


5 


0.09 ± 0.096 


5 


C/1988 Al (Liller) 


26.32 ± 0.22 


5 












-0.66 ± 0.056 


6 


0.09 ± 0.102 


5 


0.01 ± 0.071 


6 


-0.12 ± 0.043 


6 


C/2001 A2 (LINEAR) 


25.32 


1 












-0.68 


1 






0.27 


1 






C/2000 WMl (LINEAR) 


25.78 ± 0.13 


5 




2.46 


1 


0.47 ± 0.220 


3 


-0.89 ± 0.136 


5 


0.02 ± 0.169 


2 


0.02 ± 0.117 


5 






C/2()()6 VZ13 (LINEAR) 


25.35 ± 0.35 


2 


1, 


.61 ± 0.014 


2 


0.37 ± 0.093 


2 


-0.67 ± 0.028 


2 


0.15 


1 


0.10 ± 0.025 


2 






93P/Lovas 1 


24.54 ± 0.17 


2 












-0.40 ± 0.042 


2 






0.13 ± 0.375 


2 






C/2()()7 E2 (Lovejoy) 


25.07 


2 


1 


,54 ± 0.129 


2 


0.29 ± 0.093 


2 


-0.78 ± 0.035 


2 






0.07 ± 0.024 


2 






C/1994 Tl (Machholz) 


25.55 ± 0.08 


4 


2 


,21 ± 0.062 


4 


0.19 


1 


-0.62 ± 0.080 


3 






0.21 ± 0.042 


4 






C/2004 Q2 (Machholz) 


26.49 ± 0.34 


10 


1 


.96 ± 0.067 


9 


0.33 ± 0.174 


9 


-0.74 ± 0.062 


10 


0.22 ± 0.050 


10 


0.16 ± 0.035 


10 






141P/Machholz 2-A 


24.11 ± 0.06 


2 


2, 


.01 ± 0.014 


2 


0.44 ± 0.071 


2 


-0.67 ± 0.064 


2 


0.56 ± 0.144 


2 


0.13 ± 0.019 


2 






115P/Maury 


25.74 


1 












<-0.56 








<0.27 








C/1993 Yl (MR) 


25.03 ± 0.61 


7 


2 


,11 ± 0.241 


5 


0.52 ± 0.359 


7 


-0.83 ± 0.129 


7 


0.02 ± 0.414 


6 


0.05 ± 0.264 


7 






C/1980 VI (Meier) 


26.22 ± 0.47 


6 








0.25 ± 0.137 


2 


-0.71 ± 0.165 


6 


0.29 ± 0.371 


5 


0.14 ± 0.020 


5 


0.11 


1 


C/1993 Al (Mueller) 


26.18 ± 0.61 


2 




1.83 


1 


0.49 


1 


-0.35 ± 0.182 


2 


0.35 


1 


0.12 


1 






28P/Nenjmin 1 


24.62 ± 0.26 


2 












-0.17 


1 






0.27 


1 






C/1987 Bl (NTT) 


26.62 ± 0.31 


3 












-0.63 ± 0.135 


3 


0.25 ± 0.050 


2 


0.23 ± 0.176 


3 


0.11 ± 0.205 


2 


C/1992 Wl (Oshita) 


23.43 


1 












<-0.42 








0.21 


1 






C/1980 Y2 (Panther) 


26.49 ± 0.24 


4 












-0.89 ± 0.181 


3 






0.13 ± 0.077 


3 






80P/Peters-Hartley 


25.24 ± 0.14 


2 












-0.76 ± 0.122 


2 


-0.23 ± 0.057 


2 


0.15 ± 0.189 


2 


0.01 ± 0.144 


2 


83P/Russell 1 


26.13 ± 0.36 


3 












<-0.40 








<0.71 








91P/Russell 3 


24.53 


1 












<-0.04 








<0.75 








24P/Schaumasse 


24.89 ± 0.77 


9 


2, 


.48 ± 1.398 


4 


0.35 ± 0.136 


3 


-0.76 ± 0.159 


8 


0.11 ± 0.195 


3 


0.21 ± 0.349 


8 


-0.08 ± 0.272 


2 


106P/Schuster 


23.99 


1 












<-0.22 








<0.61 








31P/SW2 


25.61 ± 0.16 


6 


2, 


.17 ± 0.397 


2 


0.32 ± 0.206 


2 


-0.75 ± 0.174 


2 






-0.20 ± 0.156 


2 






102P/Shoemaker 1 


26.11 ± 0.45 


3 












-0.71 ± 0.288 


3 


0.49 


1 


0.03 ± 0.204 


3 


0.64 


1 


C/1984 Kl (Shoemaker) 


26.13 ± 0.35 


2 












-0.47 


1 






<0.51 








C/1984 U2 (Shoemaker) 


24.58 ± 0.21 


2 












<-0.84 








<-0.04 








C/1987 HI (Shoemaker) 


26.85 


1 












<0.21 








<1.00 








C/1989 A6 (Shoemaker) 


25.73 


1 












-0.31 


1 






-0.05 


1 






C/1988 Jl (SH) 


25.25 


1 












-0.57 


1 






-0.12 


1 






192P/Shoemaker-Levy 1 


24.11 


1 












<-0.72 








<0.13 








C/2002 E2 (SM) 


25.49 ± 0.04 


3 


2, 


.20 ± 0.097 


3 


0.19 ± 0.032 


3 


-0.76 ± 0.050 


3 


0.21 ± 0.021 


3 


0.20 ± 0.015 


3 






C/1986 VI (Sorrells) 


26.16 ± 0.21 


3 












-0.69 ± 0.035 


3 


0.23 ± 0.221 


2 


0.02 ± 0.062 


3 






38P/Stephan-Oterma 


26.44 ± 0.30 


14 












-0.69 ± 0.133 


12 


0.56 ± 0.231 


5 


0.03 ± 0.074 


11 






C/1983 Jl (SSF) 


24.29 ± 0.01 


2 












-0.92 ± 0.035 


2 






0.14 ± 0.057 


2 






64P/Swift-Gehrels 


24.94 ± 0.50 


10 












-0.68 ± 0.057 


10 


0.13 ± 0.064 


2 


0.14 ± 0.080 


8 


-0.07 ± 0.158 


3 


109P/Swift-Tuttle 


26.86 


1 




2.15 


1 


0.21 


1 


-0.66 


1 


0.45 


1 


0.23 


1 






C/1996 Bl (Szczepanski) 


25.91 ± 0.24 


3 


2 


,18 ± 0.053 


3 


0.31 ± 0.127 


3 


-0.75 ± 0.039 


3 


0.35 ± 0.243 


3 


0.07 ± 0.015 


3 






98P /Takamizawa 


24.54 


1 












<-0.14 








<0.65 








C/1994 32 (Takamizawa) 


26.80 ± 0.14 


2 


1 


.96 ± 0.230 


2 


0.39 ± 0.050 


2 


-0.69 ± 0.035 


2 


0.13 ± 0.272 


2 


0.08 ± 0.086 


2 






C/1994 Gl-A (TL) 


26.60 ± 0.24 


3 


2, 


.22 ± 0.046 


3 


0.33 ± 0.049 


3 


-0.69 ± 0.023 


3 


0.29 ± 0.099 


3 


0.16 ± 0.012 


3 






69P /Taylor 


24.67 ± 0.34 


5 












<-0.35 








<0.60 








9P/Tempel 1 


25.13 ± 0.21 


43 


2 


,21 ± 0.124 


31 


0.47 ± 0.144 


28 


-0.58 ± 0.334 


40 


0.31 ± 0.346 


31 


0.07 ± 0.092 


39 






lOP/Tempel 2 


25.32 ± 0.58 


26 












-0.63 ± 0.197 


18 


0.08 ± 0.100 


9 


0.05 ± 0.152 


19 


-0.03 ± 0.282 


8 


55P /Tempel-Tuttle 


25.18 ± 0.02 


2 


2, 


.14 ± 0.190 


2 


0.55 ± 0.050 


2 


-0.64 ± 0.057 


2 


0.51 ± 0.540 


2 


0.12 ± 0.014 


2 







Is3 



Comet 


log Q(CN) 

at 1 AU 


N" 








Average 


log production rate ratios 








OH/CN N" 


NH/CN 




Cs/CN 


N" 


CH/CN 


N" 


Average 
C2/CN 


N<> 


NH2/CN 




C/1987 B2 (Torasako) 


25.37 


1 








-0.78 


1 


0.02 


1 


-0.04 


1 






C1985 Tl (Thiele) 


25.69 


1 








-0.58 


1 






0.10 


1 


0.39 


1 


62P/Tsuchinshan 1 


25.10 ± 0.25 


6 








-0.70 ± 0.072 


6 


0.24 ± 0.170 


3 


0.02 ± 0.066 


6 


0.11 ± 0.067 


2 


8P/Tuttle 


25.61 ± 0.56 


2 




0.22 ± 0.246 


2 


-0.62 ± 0.122 


2 


0.33 ± 0.692 


2 


0.07 ± 0.085 


2 


-0.26 ± 0.397 


2 


41P/TGK 


24.29 


1 








-0.49 


1 






0.15 


1 






40P/Vaisala 


24.28 


1 








<-().56 








<0.27 








76P/WKI 


24.50 ± 0.06 


2 


2.19 ± 0.238 2 


0.88 ± 0.007 


2 


-0.75 


1 






-0.31 ± 0.148 


2 






81P/Wild 2 


25.55 ± 0.44 


7 








-1.05 ± 0.114 


2 






-0.27 ± 0.257 


2 






C/1986 PI (Wilson) 


26.75 ± 0.51 


4 








-0.67 ± 0.251 


3 


0.39 


1 


0.15 ± 0.150 


3 


-0.38 


1 


1 14P /Wiseman-Skiff 


24.34 ± 0.32 


3 








<-0.64 








<0.45 








43P /Wolf-Harrington 


24.97 ± 0.16 


2 








<-1.06 








<-0.08 








C/1989 Al (Yanaka) 


26.17 ± 0.57 


4 








-0.75 ± 0.236 


4 


0.42 ± 0.366 


3 


-0.10 ± 0.093 


2 







a: N = number of points used in average 

Comet names in boldface are in restricted data set; indented names are not in restricted set 

Comet Abbreviations: 



CKN — Ccrnis-Kiuchi-Nakamura 

CP — Cernis-Petrauskas 

CG — Churyumov-Gerasimenko 

dTH = Du Toit-Hartlcy 

GZ = Giacobini-Zinner 

GS = Grigg-Skjellerup 

HD = Hartley-Drinkwater 



HRC — Helin-Roman-Crockett 

HMP — Honda-Mrkos-Pajdusokova 

lAA = IRAS-Araki-Alcock 

LR = Levy-Rudenko 

Machholz 2A = Machholz Piece 2A 

MR = McNaught-Russell 



NTT — Nishikawa-Takamizawa-Tag 

SH — Shoemaker-Holt 

SSF — Sugano-Saigusa-Fujikawa 

TL = Takamizawa-Levy 

TGK = Tuttle-Giacobini-Kresak 

WKI = West-Kohoutek-Ikemura 



In Table [6] we include the comet name (some- 
times abbreviated to save space; the abbrevia- 
tions are at the bottom of the table), the aver- 
age CN production rate at 1 AU and the average 
log[Q(X)/Q(CN)] where X represents the various 
observed species. We also include the standard 
deviations of the averages (when more than one 
data point was present). The number of data 
points going into each average is listed. When 
we did not observe a species (either because it 
was not in our bandpass or we could not mea- 
sure it) the value in the table is either the upper 
limit (for C2 or C3 as described above) or is left 
blank. 

As noted above, a number of comets were ob- 
served on only one or two nights. While these 
observations may have been of extremely high 
quality (comet 122/de Vico is such an example), 
many of the single or double night observations 
were also of faint and less well observed comets. 
Thus, in a manner similar to AH95, we created 
a restricted data set for exploration of the data. 
For our purposes, we define a comet as being 
a member of the restricted data set if it was 
observed on at least three different nights (or 
night /position angle combinations). There are 
72 comets that meet this criterion. However, of 
these 72 comets, 13 had only upper limits for C2 
and/or C3. These upper limits are not constrain- 
ing and so we dropped these 13 comets from our 
definition of the restricted data set. Thus, our 
restricted data set consists of the remaining 59 
comets. In Table [6] these comets are listed with 
their names left-justified and bold-faced. The 
18 single-observation and 20 double-observation 
comets are indented and not bold-faced, as are 
the 13 comets with more than 2 observations but 
only upper limits for C3 and/or C2. This defini- 
tion of restricted data set is somewhat arbitrary 
and does not attempt to take into account the 
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Average Q(CN) at 1AU (mol sec"^) 

Figure 12: The average production rate ratios for each 
comet and each molecule are shown plotted against the 
average Q(CN) at 1 AU. In each panel, the filled symbols 
are comets in the restricted data set; the open symbols 
are the data not in the restricted data set. For C2 and 
C3, the upper limits are shown for comets with detections 
of CN but no C2 or C3 detection. The solid horizontal 
line is the weighted mean of the production rate ratios. 
The dashed lines represent the standard deviations from 
the means. For C2 and C3, the dotted lines mark the 
boundary between typical and depleted comets. Details 
of this figure are discussed in the text. 

quality of the spectra. However, it is an unbiased 
approach to choosing the restricted set. 

Figure [12] shows the production rate ratios of 
all comets as a function of the 1 AU CN produc- 
tion rate. While it would have been preferable to 
plot against CN/OH, we did not observe OH on 
enough dates to use the CN/OH data in the com- 
parison. AH95 showed that CN correlates well 
with OH (their figure 9) for all types of comets. 
The data from the restricted data set are shown 
as filled symbols, while the data that are not in 
the restricted data set are shown with open sym- 
bols The upper limits for C3 and C2 discussed 
above are shown on the appropriate panels. As 
indicated above, almost none of these limits is 
very constraining. When we did not detect a 
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species, it did not necessarily mean that there 
was none of that molecule, just that the noise 
was high. Included in Figure [12] are the weighted 
means and standard deviations for the detections 
in the restricted data set for the quantities on 
the y axes. These means and standard devia- 
tions were computed using all 59 of the comets 
in the restricted data set, without consideration 
of extreme outliers. 

The scatter in the Q(OH)/Q(CN) is quite high 
(confirmed by the standard deviation being the 
largest we have) and it is difficult to reach any 
conclusions based on these data. This is an- 
other reason why we chose not to compare all 
the comets with ratios to OH production. For 
the comets in the restricted data set, the outliers 
are 23P/Brorsen-Metcalf and C/1990 Kl (Levy) 
on the high side and 22P/Kopff and 153P/Ikeya- 
Zhang on the low side. It is Ikeya-Zhang that is 
the extremely low point. 

Q(NH)/Q(CN) is restricted to a very nar- 
row range. The three high data points are 
49P/ Arend-Rigaux, 58P/Jackson-Neujmin and 
76P/West-Kohoutek-Ikemura. These are of suf- 
ficiently good quality that they actually move 
the weighted mean upward over what the 
eye would choose. There are a number of 
comets below the bottom standard deviation 
line. C/1994 Tl (Machholz) and C/2002 E2 
(Snyder-Murakami) are the most deviant. Five 
other comets are within 0.05 of the lower stan- 
dard deviation. 

Q(NH2)/Q(CN) shows more scatter than 
Q(NH)/Q(CN), though one would assume that 
NH2 dissociates to NH. Generally, we only ob- 
served the NH2 with the IDS (with a few ex- 
ceptions) while we observed the NH with the 
LCS. Rarely are the two species observed in the 
same comet due to the wavelength range cov- 
ered in the spectra. Inspection of Table O shows 



that the error bars are generally quite a bit 
higher for the detections of NH2. The highest 
Q(NH2)/Q(CN) is for comet 102P/Shoemaker 
1, which had only one NH2 detection when the 
comet was at 2.01 AU and no NH measurement. 
27P/Crommelin and 1985 Rl (Hartley-Good) 
are slightly elevated in NH2. The lowest outliers 
are 23P/Brorsen-Metcalf, C/1980 Yl (Brad- 
field), 26P/Grigg-Skjellerup and C/1986 PI 
(Wilson). These comets had only one NH2 de- 
tection each. NH2 is a very difficult feature to 
observe since it is very weak and diffuse. Comets 
Brorsen-Metcalf and Bradfield both also had NH 
observations. Both were at the lower edge or just 
below the NH normal band. 

Q(CH)/Q(CN) seems to show many comets 
that are off the average by just about the 
standard deviation of the average. How- 
ever, there are few extreme outliers. Most 
of the extreme comets, such as C/1980 Yl 
(Bradfield), with log[Q(CH)/Q(CN)]=0.67, have 
only 1 data point. However, 2P/Encke, 
IP/Halley, C/1993 Yl (McNaught-Russeh), and 
38P/Stephan-Oterma have many data points. 
All of these comets are at the high end, except 
for McNaught-Russell. 23P /Brorsen-Metcalf has 
three data points and is also high. The low- 
est data point is comet 80P/Peters-Hartley, with 
two data points when the comet was at 1.78 AU. 

If we ignore the upper limits, Q(C2)/Q(CN) 
shows a very tight correlation for most comets 
but some comets show significant devia- 
tions. The high data point is 68P/Klemola. 
The low outliers are comets 19P/Borrelly, 
32P/Comas Sola, 21P/Giacobini-Zinner, 
75P/Kohoutek, C/1988 Jl (Shoemaker-Holt), 
31P/Schwassmann-Wachmann 2, 76P/West- 
Kohoutek-Ikemura, 81P/Wild 2, and C/1989 Al 
(Yanaka). Most of these are well observed. The 
two lowest comets are Giacobini-Zinner (-0.57) 
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and Comas Sola (-0.89). 

Q(C3)/Q(CN) also shows high scatter, though 
most of the upper limits are not very con- 
straining. The high outliers are 49P/Arend- 
Rigaux, C/1995 01 (Hale-Bopp), 88P/Howell, 
58P/Jackson-Neujmin, C/1990 Kl (Levy), and 
28P/Neujmin 1. Neujmin 1, Jackson-Neujmin 
and Howell are the three highest values. On 
the low side, 21P/Giacobini-Zinner is the ex- 
tremely low point, with log[Q(C3)/Q(CN)]=- 
1.47. The other low points (values less than 
-0.9) are 19P/Borrelly, 153P/Ikeya-Zhang, and 
81P/Wild 2. 

Categorizing individual comets was not the 
point of our program. We sought to understand 
whether comets could have been formed with 
more than one compositional type. This is not 
the first study of its kind. In addition to the 
papers by AH95, F09 and LSll listed in the in- 
troduction, smaller studies were carried out by 
A'Hearn and Millis (1980), Newburn and Spin- 
rad (1984), Fink and Hicks (1996) and by our 
group (Cochran 1987; Cochran et al. 1987, 1989, 
1992; Cochran and Barker 1987). Thus, it has 
already been discovered that, while the spectra 
of most comets are very similar, some comets 
show much weaker C2 and C3 relative to other 
species. AH95 used a database of 85 comets to 
quantify this trend and have defined two taxo- 
nomic classes: typical and carbon-chain depleted 
(shortened to depleted). 

In order to study the abundance correlations 
of the various species, we started by forming 
weighted means and standard deviations of the 
data. These are listed in Table [71 In this ta- 
ble, we list the values for each species over var- 
ious subsets of the data. First, we computed 
the mean for all 110 comets with CN observa- 
tions (upper limits were not used in the means). 
We also list the means and standard deviations 



for the complete restricted data set. Inspection 
of Table [7] shows that the restricted data set is 
representative of the whole data set. The com- 
plete restricted data set's means and standard 
deviations are shown as solid and dashed lines, 
respectively, in Figure [T2l 

Using the values from the complete data set, 
we have formed a definition for a depleted 
comet in our data set. A depleted comet is 
one for which log[Q{C^)/Q{CN)] < -0.86 and 
log[Q(C2)/Q(CiV)] < 0.02. Note that we re- 
quire both the C3 and C2 ratios to be low. 
These values are shown in Figure [12] as dot- 
ted lines in the C2 and C3 panels. Within the 
restricted data set we find 5 comets that fit 
this definition: 19P/Borrelly, 32P/Comas Sola, 
21P/Giacobini-Zinner, C/2000 WMl (LIN- 
EAR), and 81P/Wild 2. Borrelly, Giacobini- 
Zinner and Wild 2 have been previously shown 
to be depleted by us and others. Thus, in our 
restricted data set we find 5/59 comets belong 
to the carbon-chain depleted group, or 9%. For 
the comets not in the restricted set, C/1979 Yl 
(Bradfield) is also depleted. In addition, while 
most of the upper limits are not constrain- 
ing, two of the comets with only upper limits 
have values that are below the defined cutoff 
for depleted comets and therefore must be de- 
pleted: C/1984 U2 (Shoemaker) and 43P/Wolf- 
Harrington (both were found depleted by AH95). 

For these depleted comets, the depletion of C2 
and C3 is not a subtle effect. Figure [13] shows 
spectra of 21P/Giacobini-Zinner and 8P/Tuttle, 
obtained with the same instrument. The view- 
ing geometry for these two observations is almost 
identical, with only a very slight difference in he- 
liocentric distance. The two spectra are scaled 
by the CN band. The Tuttle spectrum shows 
all of the usual molecules while the Giacobini- 
Zinner spectrum is "missing" the C2 and C3. In- 
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Table 7: Derived Production Rate Ratios 



Log Production Rate Ratio with respect to CN 





OH 


NH 


Ca 


CH 


Ca 


NHa 












(average) 




All Comets 


2.09 ± 0.35 


0.46 ±0.19 


-0.68 ±0.20 


0.25 ±0.19 


0.15 ±0.20 


-0.09 ±0.27 


(110 comets) 


(33 comets) 


(38 comets) 


(84 comets) 


(41 comets) 


(84 comets) 


(29 comets) 


All Restricted Comets 


2.29 ±0.32 


0.48 ±0.21 


-0.68 ±0.19 


0.25 ±0.18 


0.15 ±0.23 


-0.09 ±0.25 


(59 comets) 


(23 comets) 


(27 comets) 


(59 comets) 


(41 comets) 


(59 comets) 


(26 comets) 


Restricted, Typical 


2.29 ±0.33 


0.48 ±0.21 


-0.67 ±0.15 


0.25 ±0.18 


0.16 ±0.14 


-0.09 ±0.26 


(54 comets) 


(21 comets) 


(25 comets) 


(54 comets) 


(39 comets) 


(54 comets) 


(25 comets) 


Restricted, Depleted 


2.19 ±0.16 


0.44 ± 0.28 


-1.06 ±0.39 


0.02 ± 0.26 


-0.24 ±0.60 


0.01 


(5 comets) 


(2 comets) 


(2 comets) 


(5 comets) 


(2 comets) 


(5 comets) 


(1 comet) 



Depleted comet is defined as log [Q(C3)/Q(CN)] < -0.86 and log [Q(C2)/Q(CN)] < 0.02. 



spection of the Giacobini-Zinner spectrum shows 
we detected the weak CN Af = —1 band at 
~ 4200 A, so the failure to detect the C2 and 
C3 is not a failure of signal/noise. Giacobini- 
Zinner really does have very little of these two 



21 P/Giacobini-Zinner 
23 Sept 1 985 
A=0.51au R^=1.06au 




8P/Tuttle 
9 Dec 1980 
A=0.51au R^=1.02au 



-I 
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Wavelength (A) 

Figure 13: IDS Spectra of 21P/Giacobini-Zinner and 
8P/Tuttle are shown when both comets were near a he- 
liocentric distance of 1 AU and a geocentric distance of 
0.5 AU. While both comets show similar CN (including 
the weaker A« = — 1 band near ~4200A), there is little 
evidence of the C3 or C2 bands in the Giacobini-Zinner 
spectrum. 



molecules! 

In Table [7| we list means and standard devia- 
tions for the 54 comets in the restricted data set 
that are not depleted ( "typical" ) and separately 
list these values for the 5 depleted comets. Com- 
parison of these two groups shows that, while 
the depleted comets are significantly different in 
their C3 and C2 ratios, OH and NH look the 
same in both sets. CH shows a little difference 
between the typical and depleted comets, but 
the error bars are high. NH2 has too little data 
(only 1 comet) in the depleted group to draw any 
conclusion. 

As noted above, other groups have looked at 
abundance correlations between comets. As each 
study has used slightly different fluorescence effi- 
ciencies and Haser scale lengths, one cannot sim- 
ply compare the production rate ratios from one 
study with another study. However, as AH95 
pointed out "...the pattern of values of the ra- 
tios will, to first order, be independent of these 
parameters" . 

Figure [T3] illustrates the effects of using differ- 
ent scale lengths for computing the column den- 
sities. It shows the CN and C2 column densities 
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for two different comets along with model fits 
using the scale lengths of AH95, F09, LSll and 
the values used in this paper. These two particu- 
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log Position (km) 

Figure 14: The CN and C2 data for comets C/1988 Al 
(Liller) and C/1996 Bl (Szczepanski) are shown along 
with Haser model fits using the scale lengths of this pa- 
per, AH95, F09 and LSll. The two comets were observed 
at comparable heliocentric distances. For C2, the fits us- 
ing our scale lengths are comparable to those using the 
scale lengths of F09. The AH95 and LSll fits are similar. 
For CN, our fits are indistinguishable from ones using the 
scale lengths of LSll. The AH95 and F09 models are 
higher and lower, respectively. However, note the differ- 
ent production rates that are computed with these fits. 
See text for discussion. 

lar cometary data sets were chosen because they 
are LCS data with lots of points, they are high 
quality so that the trends are well defined, there 



is very little asymmetry seen and they are both 
at approximately the same heliocentric distance. 
This distance is sufficiently far from 1 AU that 
the various heliocentric scaling laws play a sig- 
nificant role in the scale lengths. C2 was chosen 
because it is in the bandpass of all of the surveys. 
CN is the molecule used in the denominator of 
the ratios. Inspection of the two C2 panels shows 
that our scale lengths do a reasonably good job 
of fitting the comet Liller data, though they may 
underestimate the column densities around 10^ 
km from the optocenter. For comet Szczepan- 
ski, the fit is good at cometocentric distances 
greater than 1000 km but the inner part of the 
coma is fit poorly. This illustrates that the Haser 
model does not handle the two-step (grandpar- 
ent to parent to daughter) dissociation well. The 
models using the F09 scale lengths are almost in- 
distinguishable in shape to those with our scale 
lengths, but the derived log production rates 
differ from ours by ~0.25. The model fits us- 
ing the AH95 and LSll scale lengths are very 
similar. However, the production rates for the 
AH95 scale lengths are essentially the same as 
ours while the LSll production rates are inter- 
mediate to ours and F09. For CN, our model fits 
and those of LSll arc indistinguishable but yield 
production rates that differ by 0.12. Models us- 
ing F09 and AH95 scale lengths are below and 
above our model fits, respectively. Again, our 
production rates match those derived from AH95 
scale lengths, while LSll and F09 scale lengths 
yield the same, higher, production rates. For 
Szczepanski, our CN scale lengths fit extremely 
well, while we overestimate the column densities 
in the inner coma of Liller. 

The scale lengths used for each model are 
on each panel and are considerably different for 
the different groups. It is a well known ef- 
fect that there is a family of Haser model scale 
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lengths allowed for each molecule since the par- 
ent and daughter scale lengths are degenerate in 
the model and can compensate for one another. 
There is nothing in Figure O to warrant picking 
another set of scale lengths over the ones we have 
used in this paper (Table [2]) . What is a bit more 
critical is the boundary for declaring a comet to 
be typical or depleted. We would say that Liller 
was slightly depleted. AH95 listed Liller, though 
not in their restricted set, as typical. Using their 
scale lengths and our cutoff, we would still find 
Liller to be depleted. With the AH95 definition 
and either our or their scale lengths Liller would 
be typical. We found Szczepanski to be typical, 
as did F09. 

The recent study of LSll has far fewer comets 
than our study or the others we compare to. 
However, the data should be comparable to ours 
as they used a long slit CCD instrument with 
a bandpass similar to the LCS. They suffered 
from difficult sky subtraction because of street 
lamps near Lick Observatory, something they 
took great pains to deal with. McDonald does 
not have this problem. Also, some of their ob- 
servations were at twilight when the sky contri- 
bution was more severe. We almost never had 
to observe under such conditions so did not en- 
counter these problems. 

F09 used a long slit CCD instrument but con- 
centrated on the red region of the spectrum. 
Thus, he was unable to detect OH, NH or C3. 
For CN, he observed the CN red system instead 
of the violet system we observed. F09 observed 
NH2 and O (^D). The latter could be used to 
derive QlHgO). 

AH95 used a fundamentally different data 
type. Their survey used photometry and narrow 
band filters optimized for comets. With photom- 
etry, the apertures are considerably larger than 
with our spectrograph and so AH95 are much 



more sensitive at low signal levels. They cannot, 
however, detect the very weak bands, such as 
CH. In addition, since the different species are 
observed sequentially and not simultaneously, 
they require photometric conditions. Analyses 
of photometric observations must assume spher- 
ical symmetry of the coma since they get little 
spatial information. 

AH95 looked at 85 comets to determine abun- 
dance patterns; 41 of those comets formed their 
restricted data set over which they drew con- 
clusions. They found that 12/41 (29%) of the 
comets were in the depleted group. This is sig- 
nificantly more than we found. However, AH96 
used a definition based only on a low value for 
C2. Some of their depleted comets have low C2 
but show C3 within normal bounds. This can be 
seen clearly in their Figure 10, comparing panels 
a and b. 

Using only a criterion for depleted comets 
that log[Q{C2)/Q{CN)] < 0.02, with no limit 
on the C3, we add ten comets to our list 
of depleted comets: 71P/Clark, 26P/Grigg- 
Skjellerup, 75P/Kohoutek, C/1988 Al (Liller), 
31P/Schwassmann-Wachmann 2, C/1988 Jl 
(Shoemaker-Hoh), C/1986 VI (Sorrells), 
62P/Tsuchinshan 1, 76P/ West-Kohoutek- 
Ikemura, and C/1989 Al (Yanaka). With this 
definition, we find 15/59 depleted comets, or 
25%. This is consistent with the quantity found 
by AH95. These 10 comets all show ratios for 
C3 production at an average value or below. 
However, none are close to our defined value for 
C3 depletion. 

AH95 also noted that some comets were en- 
hanced in C2/CN. We too see such a trend, al- 
though not necessarily the same comets. We 
both found C/1980 Yl (Bradfield) to be en- 
hanced in C2. In addition, 68P/Klemola has a 
strong C2 enhancement. 
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F09 did not observe C3 in their bandpass 
and so their definition of depleted comets was 
based entirely on the C2 abundance. F09 ob- 
served a total of 92 comets, with 50 in their 
restricted group. F09 broke their comets into 
more groups than just typical and depleted. 
However, all comets not in their typical group 
had C2 depleted relative to CN (except possibly 
C/1988 Yl (Yanaka)). They found 30% depleted 
comets. 

LSll observed 26 comets. They did not define 
a subset of comets as part of a restricted group. 
Many of their comets were observed on only 1 
night. LSll chose to derive new Haser model 
scale lengths for their data set because they did 
not feel that any of the published scale lengths fit 
their data well. In addition, they derived a he- 
liocentric distance scaling for those scale lengths 
that was different (not RfJ than most researchers 
use. It was not clear if they used mean scale 
lengths or used different scale lengths for each 
comet. In addition, it appears from their fig- 
ures that they did not remove the underlying 
solar continuum before computing the column 
densities, though they did fit a continuum be- 
fore integration. For dusty comets, this will in- 
crease the difficulty of detecting the gas. Once 
they had computed their production rates in this 
manner, LSll used the AH95 definition of de- 
pleted comets, finding all their Jupiter family 
comets depleted (including 4P/Faye where they 
did not detect any C2 and 46P/Wirtanen that 
has a much higher C2/CN ratio in their data 
than this cutoff). They also indicated that as 
a class the long period comets were borderline 
depleted. LSll did not attempt to convert the 
AII95 limits to their scale lengths so it is not 
obvious that the comets they claim to be de- 
pleted would be depleted if they derived a cut- 
off from their data using their scale lengths. Us- 



ing the AH95 definition and the data from their 
Table 11, there are 13/26 comets that are de- 
pleted, or 50% (this includes several with no C2 
detection that they take to indicate a comet is 
depleted). 

Table [8] lists the 110 comets we observed that 
had at least CN detected, along with a sum- 
mary of our findings. As before, the restricted 
comets are listed in bold and the remainder are 
indented. The group of comets that are depleted 
under our most restrictive definition of low C2 
and C3 are noted as "Depleted" (with the word 
depleted underlined). The additional 10 comets 
with low C2 but not low C3 are listed as "(de- 
pleted)". Comets that have all of their produc- 
tion rate ratios within the normal ranges shown 
for the restricted set in Table [7] are marked "Typ- 
ical" . There are comets that are not depleted but 
have at least one production rate ratio outside of 
the normal ranges. We note the species that are 
either high or low. Inspection of the table shows 
that some comets have only one species outside 
a normal range while others (e.g. C/1990 Yl 
(Bradfield)) have many species outside the nor- 
mal range. Comets with only limits are so noted. 

Included in Table [8] are the findings of the 
other studies discussed above for comets ob- 
served in common. When comparing our comets 
with the findings of others, it is important to re- 
member that only C2 was used in other studies as 
a criterion for being called depleted, so that any 
comet where we saw normal C2 but other species 
out of range would still be typical in other stud- 
ies. 

Inspection of Table [8] shows that we are gen- 
erally in (remarkably good) agreement on how 
we classify comets, with some exceptions. As 
noted above when we examined scale lengths, 
placement of the cut-off for defining depletion 
has an important affect on our classifications. 
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Table 8: Summary of Abundance Patterns 



Comet 


Dynamical 


This 


A'Hearn 


Fink 


Langlund-Shula 




Type 


Paper 


et ai. 




& Smith 








(1995) 


(2009) 


(2011) 


49P / Arend-Rigaux 


JFC 


high NH, C3 


Typical 






C/1982 Ml (Austin) 


LPC 


Typical 


Enhanced 






C/1984 Nl (Austin) 


LPC 


Typical 


Typical 






C/1989 XI (Austin) 


LPC 


low NH 


Typical 


Typical 




19P/Borrelly 


JFC 


Depleted + low CH 


Depleted 


low C2 




140P/Bowcll-Skifr 


JFC 


low C3, high CH 








C/1979 Yl (Bradficld) 


LPC 


Depicted 


Typical 






C/1980 Yl (Bradfleld) 


LPC 


high CH, C2; low NH2, C3 


Enhanced 






C/1987 PI (Bradfleld) 


LPC 


Typical 


Typical 






16P/Brooks 2 


JFC 


Limits 


Depleted 






23P /Brorsen-Metcalf 


HTC 


high OH, CH; low NH, NH2 


Typical 


Typical 




C/1983 Ol (Ccrnis) 


LPC 


high Cs 


Depleted 






C/1990 El (CKN) 


LPC 


low CH 








C/1980 Ol (CP) 


LPC 


Limits 


Typical 






67P/CG 


JFC 


Typical 


Depleted 


low C2 




71P/Clark 


JFC 


low C2 (Depleted) 


low C2 






32P/Comas Sola 


JFC 


Depleted 








27P/Crommelin 


HTC 


low NH2 


Typical 






6P/d' Arrest 


JFC 


Typical 


Typical 


Typical 




122P/dc Vico 


HTC 


low C3 


Typical 






79P/dTH 


JFC 


Limits 








2P/Encke 


JFC 


low NH, high CH 


Typical 


Typical 




4P/Faye 


JFC 


Limits 


Depleted 




Depleted 


37P/Forbes 


JFC 


Limits 








78P/Gehrels 2 


JFC 


Limits 


Typical 






21P/GZ 


JFC 


Depleted 


Depleted 


low C2, NH2 




84P/Giclas 


JFC 


Limits 








26P/GS 


JFC 


low C2, NH2 (Depleted) 


Depleted 






65P/Gunn 


JFC 


Limits 


Depleted 






C/1995 Ol (Hale-Bopp) 


LPC 


high C3, CH 




Typical 


Depleted 


IP/Halley 


HTC 


high CH, C2 


Typical 


Typical 




51P /Harrington 


JFC 


Typical 








52P /Harrington- Abell 


JFC 


Limits 








C/1995 Q2 (HD) 


LPC 


Limits 








1985 Rl (Hart ley- Good) 


LPC 


high NH2 


Typical 


Typical 




161P/Hartlcy-IRAS 


JFC 


high NH2 


Typical 






lllP/HRC 


JFC 


Limits 








45P/HMP 


JFC 


low C3 


Typical 


low C2 




88P/Howell 


JFC 


high C3 


Typical 






C/1995 Yl (Hyakutake) 


LPC 


low OH, NH 




Typical 




C/1996 B2 (Hyakutake) 


LPC 


Typical 




Typical 




153P/Ikeya-Zhang 


LPC 


low OH, C3; high CH 




Typical 




C/1983 HI (lAA) 


HTC 


Typical 


Typical 






126P/IRAS 


JFC 


Limits 


Depleted 






C/1983 02 (IRAS) 


LPC 


Limits 


Typical 






58P / Jackson- Neujmin 


JFC 


high NH, C3 




Typical 




59P /Kearns-Kwee 


JFC 


Typical 


Typical 






68P/Klemola 


JFC 


high C2 


Depleted 
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Table [8] continued 



Comet 


Dynamical 


This 


A'Hearn 


Fink 


Langlund-Shula 




Type 


Paper 


et al.. 




& Smith 








(1995) 


(2009) 


(2011) 


75P/Kohoutek 


JFC 


low C2, NH2 (Depleted) 








22P/Kopff 


JFC 


low OH 


Typical 


Typical 


Depleted 


144P/Kushida 


JFC 


Typical 








C/1987 Al (Levy) 


LPC 


low CH 








C/1990 Kl (Levy) 


LPC 


high OH, C3, C2 


Typical 


Typical 




C/1984 VI (LR) 


LPC 


Typical 


Typical 






C/1988 Al (Liller) 


LPC 


low C2 (Depleted) 


Typical 






C/2001 A2 (LINEAR) 


LPC 


Typical 








C/2000 WMl (LINEAR) 


LPC 


Depleted + low CH 








C/2006 VZ13 (LINEAR) 


LPC 


low OH 








93P/Lovas 1 


JFC 


high C3 




Typical 


Depleted 


C/2007 E2 (Lovcjoy) 


LPC 


low OH 






Typical 


C/1994 Tl (Machholz) 


LPC 


low NH 




low C2 




C/2004 Q2 (Machholz) 


LPC 


Typical 








141P/Machholz 2-A 


JFC 


high CH 




Typical 




115P/Maury 


JFC 


Limits 








C/1993 Yl (MR) 


LPC 


low CH 




Typical 




C/1980 VI (Meier) 


LPC 


low NH 


Typical 






C/1993 Al (Mueller) 


LPC 


low OH, high C3 








28P/Neujmin 1 


JFC 


high C3 


Typical 






C/1987 Bl (NTT) 


LPC 


Typical 


Typical 






C/1992 Wl (Oshita) 


LPC 


limits 








C/1980 Y2 (Panther) 


LPC 


low C3 


Typical 






80P/Peters-Hartley 


JFC 


low CH 




Typical 




83P/Russell 1 


JFC 


limits 








91P/Russell 3 


JFC 


limits 








24P/Schaumasse 


JFC 


Typical 




Typical 




106P/Schuster 


JFC 


limits 








31P/SW2 


JFC 


low C2 (Depleted) 


Depleted 






102P/Shoemaker 1 


JFC 


high CH, NH2 








C/1984 Kl (Shoemaker) 


LPC 


high C3 


Enhanced 






C/1984 U2 (Shoemaker) 


LPC 


limits/Depleted 


Depleted 






C/1987 HI (Shoemaker) 


LPC 


limits 








C/1989 A6 (Shoemaker) 


LPC 


high C3, low C2 (Depleted) 








C/1988 Jl (SH) 


LPC 


low C2 (Depleted) 








192P/Shoemaker-Lcvy 1 


JFC 


limits 








C/2002 E2 (SM) 


LPC 


low NH 








C/1986 VI (Sorrells) 


LPC 


low C2, NH2 (Depleted) 


Depleted 






38P/Stephan-Oterma 


HTC 


high CH 


Typical 






C/1983 Jl (SSF) 


LPC 


low C3 


Typical 






64P/Swift-Gehrels 


JFC 


Typical 








109P/Swift-Tuttle 


HTC 


low NH, CH 


Typical 


Typical 




C/1996 Bl (Szczepanski) 


LPC 


Typical 




Typical 




98P/Takamizawa 


JFC 


Limits 


Depicted 






C/1994 J2 (Takamizawa) 


LPC 


Typical 




low C2 




C/1994 Gl-A (TL) 


LPC 


Typical 








69P/Taylor 


JFC 


limits 








9P/Tempel 1 


JFC 


Typical 


Typical 


low C2 


Depleted 



34 



Table |8] continued 



Comet 


Dynamical 


This 


A'Hearn 


Fink 


Langlund-Shula 




type 


Paper 


a+ ft] 




cz omitn 








(1995) 


(2009) 


(2011) 


lOP/Tempel 2 


JFC 


Typical 


Typical 


Typical 


Depleted 


55P/Tempel-Tuttle 


HTC 


high CH 




Typical 




C/1987 B2 (Terasako) 


LPC 


low CH, C2 (Depleted) 








C1985 Tl (Thielc) 


LPC 


high NH2 


Typical 






62P/Tsuchinshan 1 


JFC 


low C2 (Depleted) 


Typical 


Typical 




8P/Tuttle 


HTC 


low NH 


Typical 




Typical 


41P/TGK 


JFC 


high C3 








40P/Vaisala 


JFC 


limits 








76P/WKI 


JFC 


high NH, low C2 (Depleted) 








81P/Wild 2 


JFC 


Depleted 


Depleted 


low C2 


Depleted 


C/1986 PI (Wilson) 


LPC 


low NH2 








114P/Wiseman-Skifr 


JFC 


limits 








43P/Wolf-Harrington 


JFC 


Limits /Depleted 


Depleted 


low C2, NH2 




C/1989 Al (Yanaka) 


LPC 


low C2 (Depleted) 




Typical 





In the discussion here, keep in mind our reser- 
vations about the definition of depletion used 
by LSI 1. In our restricted data set, we find 
C/1982 Ml (Austin) typical while AH95 find it 
enhanced. F09 and AII95 find 67P/Churyumov- 
Gerasimenko depleted but we find it typical. 
C/1995 01 (Hale-Bopp) is typical in our data 
and that of F09 but found to be depleted 
in LSll. AH95 find 68P/Klemola to be de- 
pleted but we find it to be enhanced in C2. 
22P/Kopff and lOP/Tempel 2 are typical for us, 
F09 and AH95, but LSll find them depleted. 
9P/Tempel 1 appears normal to us and AII95, 
but depleted to F09 and LSll. C/1988 Al 
(Liller), 62P/Tsuchinshan 1, and C/1989 Al 
(Yanaka) are depleted in our extended definition 
(Liller and Tsuchinshan 1 barely so) but typi- 
cal according to AH95 and/or F09. C/1994 Tl 
(Machholz) is depleted for F09 but typical for 
us. In addition, there are five comets not in our 
restricted set where we see some differences with 
others. 

If all of the depleted comets started in the 
same formation reservoir, we would expect to 
see all of the depleted comets having the same 



dynamical type. Of course, the dynamical type 
does not necessarily correspond to a single reser- 
voir since there have been exchanges of objects 
from one reservoir to another via gravitational 
perturbations. For our 59 comets in the re- 
stricted data set, 27 are JFCs, 5 are HTCs and 27 
are LPCs. Using our strict definition of carbon- 
chain depleted comets, we find 4/5 of them are 
JFCs (80%) and 1 is an LPC (20%). Using the 
more relaxed definition requiring only C2 to be 
depleted, we find 10 of the depleted comets are 
JFCs (67%) and 5 are LPCs (33%). For a partic- 
ular dynamical type and the relaxed definition, 
10/27 or 37% of the JFCs are depleted; 5/27 
or 18.5% of the LPCs are depleted. We see no 
depleted HTCs, but with only 5 HTCs in our 
sample this might not be significant. 

In Table O we show the average production 
rate ratios function of dynamical type, com- 
bining depleted and typcial comets in each dy- 
namical type. Inspection of Table [9] shows that 
the averages are indistinguishable for OH (with 
only 1 Halley Type comet with OH the devia- 
tion of this group is not meaningful), or C3. CH 
and NH2 might show differences with dynami- 
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Table 9; Production Rate Ratios as a function of Dynamical Type 



Log Production Rate Ratio with respect to CN 





OH 


NH 


Ca 


CH 


Ca 
(average) 




Jupiter Family 


2.15 ±0.32 


0.60 ±0.24 


-0.69 ± 0.25 


0.11 ±0.19 


0.04 ±0.31 


-0.03 ±0.28 


(27 comets) 


(11 comets) 


(12 comets) 


(27 comets) 


(14 comets) 


(27 comets) 


(11 comets) 


Halley Type 


2.76 


0.25 


-0.67 ±0.08 


0.46 ±0.28 


0.16 ±0.13 


0.10 ±0.34 


(5 comets) 


(1 comet) 


(1 comet) 


(5 comets) 


(4 comets) 


(5 comets) 


(4 comets) 


Long Period 


2.17 ±0.38 


0.36 ±0.18 


-0.68 ±0.12 


0.27 ±0.17 


0.18 ±0.11 


-0.17 ±0.22 


(27 comets) 


(11 comets) 


(14 comets) 


(27 comets) 


(23 comets) 


(27 comets) 


(11 comets) 



cal type but the error bars are large (and over- 
lapping) and there is no obvious sequence with 
current orbit size. For NH, there appears to be 
a significant difference between JFCs and LPCs 
(again HTCs do not have enough data to be 
meaningful) but the error bars are overlapping. 
For C2, there is a significant difference between 
JFCs and other dynamical types, though LPCs 
and HTCs are identical. This trend for C2 is 
consistent with us finding many more JFCs that 
are depleted and is in agreement with the trend 
found by AH95. We concur with their conclusion 
that the compositional similarities between dy- 
namical types are suggestive that the interiors of 
comets are similar to their exteriors. In contrast, 
LSll reported a trend that the C2, NH and NH2 
increased in a regular fashion as the comet or- 
bits became larger (though their Table 12 shows 
the C2/CN decreasing from HTCs to LPCs and 
being the same for dynamically new and HTCs) . 

Perhaps the differences we see for a large 
number of JFCs is the result of evolution- 
ary differences instead of formational differ- 
ences. While we cannot categorically refute 
this concept, we point to observations of comet 
73P/Schwassmann-Wachmann 3 (SW3) as evi- 
dence to the contrary. Comet SW3 split into 
several pieces in 1995. It made a close ap- 



proach to the Earth in 2003 and was exten- 
sively studied. 73P/Schwassmann-Wachmann 3 
showed no compositional differences between dif- 
ferent pieces (Kobayashi et al. 2007, Jehin et al. 
2008, Weaver et al. 2008, Schleicher and Bair 
2011), suggesting that the comet was the same 
composition throughout. SW3 is a JFC that is 
depleted. Therefore, its depletion is not just a 
surface effect. 

AH95 found a trend of dust-to-gas ratio as a 
function of cometary perihelion distance (their 
Figure 4). With our small apertures, we do not 
feel comfortable computing A//? so cannot com- 
ment on whether our data show such a trend. 
We did, however, examine the production rate 
ratios function of cometary perihelion dis- 
tance. These are shown in Figure [T5j We find 
no believable trend of the production rate ratios 
with perihelion distance. 

5. Summary 

When we began this survey in 1980, we fully 
expected that we would find some comets with 
very different spectra than the vast majority of 
comets. Instead, we found a mostly homoge- 
neous group of objects with subtle compositional 
differences. In particular, we found: 
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Figure 15: The production rate ratios are plotted as 
a function of the periheUon distances of the current 
cometary orbits. The symbols used are the same as in 
Figure 1121 as are the definitions of the horizontal lines. 
There are no believable trends of the ratios as a function 
of perihelion distance. 



• The vast majority of comets have remark- 
ably similar compositions. 

• There exists a group of comets which are de- 
pleted in the carbon-chain molecules. Using 
a very strict definition of carbon-chain de- 
pletion requiring depletion in both C2 and 
C3, we find that 9% of the comets are de- 
pleted. Using a more relaxed definition that 
requires only the C2 be depleted (similar to 
other authors), we find 25% of the comets 
are depleted. This is consistent with other 
authors. 

• For the most part, we agree with other stud- 
ies on which individual comets show deple- 
tions. 

• Carbon-chain depleted comets can be of any 
dynamical type (though we did not observe 
any HTCs with depletion our HTC sample 
is very small). 



Using the relaxed definition, two thirds of 
the depleted comets are JFCs and one third 
are LPCs (with the stricter definition it is 
80% and 20% respectively). 

Depleted comets make up 37% of the JFCs 
and 18.5% of the LPCs. 

On average, there are few differences in com- 
position with increasing current orbit semi- 
major axis. However, C2 and NH may show 
significant trends of lower C2 and higher NH 
for JFCs than for LPCs. 
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